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ABSTRACT 


This  report  describes  the  procedure  utilized  in  the  assembly  and 
evaluation  of  a  small  intermittent  blowdown  supersonic  wind  tunnel. 

The  tunnel  proper  was  fully  constructed  and  nozzle  blocks  for 
attainment  of  a  Mach  number  of  1.5  or  2.37  in  the  1/2”  by  1.718” 
test  section  were  available.  The  basic  schlieren  system  was  purchased 
complete.  An  attempt  was  made  to  use  or  convert  readily  available 
materials  for  all  components. 

It  was  expected  the  results  of  this  project  would  be  of  a  perman¬ 
ent  nature  and  would  be  utilized  by  the  Mechanical  Engineering  Depart - 

•% 

ment  as  a  laboratory  installation.  The  installation  is  suitable  for 
the  study  of  compressible  fluid  flow  and  the  phenonemon  of  normal 
shocks  but  is  not  suitable  for  the  study  of  the  performance  of  shapes 
or  models  in  the  fluid  flow. 

The  tunnel  installation  operates  satisfactorily  with  a  maximum 
Mach  number  of  2.29  with  approximately  0.28  lbm  of  air  per  second. 
Schlieren  photographs  of  the  flow  field  may  be  taken  in  either  black 
and  white  or  in  color. 


ii 


ACKNOWLEDGEMENT 


The  author  wishes  to  express  his  appreciation  to  Professor  Charles 
P.  Howard  for  his  guidance,  assistance  and  numerous  benefical  sugges¬ 
tions  during  all  phases  of  the  project.  Appreciation  is  also  extended 
to  Messrs.  K.  W.  Mothersell  and  Joseph  Beck  for  their  assistance  in 
the  construction  and  installation  of  the  tunnel  components. 


LIST  OF  ILLUSTRATIONS 

Figure  Title  Page 

1  Schematic  Diagram  of  the  Tunnel  System  22 

2  Test  Section  Subassembly  23 

3  Compressor  Discharge  Pressure  Control  and 

Air  Drying  System  24 

4  Tunnel  Section  Stagnation  Pressure  Control 

System  24 

5  Stilling  Tank  25 

6  Test  Section,  Disassembled  26 

7  Stilling  Tank,  Entrance  Nozzle,  Test  Section 

and  Diffuser  Assembly  27 

8  Overall  View  of ’Tunnel  Section  from  the  Rear 

Side  29 

9  Overall  View  of  Tunnel  Section  from  the 

Front  Side  28 

10  Manometer  Manifolds  25 

11  Schematic  Diagram  of  Schlieren  System  30 

12  Schlieren  System  Light  Source,  Knife  Edge, 

and  Camera  Viewing  Components  31 

13  Schlieren  Photographs  of  Tunnel  Running 

Shock  Free  32 

14  Schlieren  Photographs  Demonstrating  Effect 
of  Increased  Back  Pressure  on  the  Normal 

Shock  Location  33 

15  Schlieren  Photographs  Demonstrating  Effect 

of  Reducing  the  Stagnation  Pressure  on  the 
Location  of  a  Normal  Shock  in  the  Nozzle  34 

16  Schlieren  Photograph  Demonstrating  Effect 

of  Minimum  Starting  Area  of  Tunnel  35 

17  Schlieren  Photographs  of  a  Normal  Shock  in 
the  Diffuser  Section  and  Effect  of  Oil 

Vapor  in  the  Air  36 

18  Schlieren  Photographs  Demonstrating  Mach 

Number  Determination  37 

iv 


Figure  Title  Page 

19a  Nozzle  Block  Contours  vs.  Nozzle  Length  38 

19b  Diffuser  Contours  vs.  Nozzle  Length  38 

20  Adiabatic  Mach  Number  vs.  Tunnel  Length 

with  Tunnel  Running  Shock  Free  and  Dif¬ 
fuser  Set  at  Constant  1.718  Inch  Height  38 

21  Adiabatic  Mach  Number  vs.  Tunnel  Length 

with  Tunnel  Running  Shock  Free  and  Dif¬ 
fuser  Set  at  Maximum  1.883  Inch  Height  38 

22  Adiabatic  Mach  Number  vs.  Tunnel  Length 

with  a  Normal  Shock  Standing  in  Nozzle  39 

23  Tunnel  Pressure  Ratios  vs.  Tunnel  Length 
with  and  without  a  Normal  Shock  Standing 

in  Nozzle  39 

24  Adiabatic  Mach  dumber  vs.  Tunnel  Length 
Tunnel  Operating  (1)  with  Increased  Back 
Pressure  and  (2)  Diffuser  Set  at  Minimum 

1.40  Inches  39 

25  Orifice  Discharge  Coefficient  vs.  Stilling 

Tank  Stagnation  Pressure  39 

1-A  Pressure  Differential  Between  Stagnation 
Pressure  and  Tunnel  Test  Section  Pressure 
Required  to  Operate  Tunnel  Supersoni¬ 
cally  47 

1-C  Flow  Coefficient  K  vs.  Reynolds  Number 

for  Flat  Plate  Orifice  53 

1- D  Stress  Distribution  of  Diaphragm  with 

Clamped  Edges  56 

2- D  Exterior  of  Transducer  58 

3- D  Interior  of  Transducer  58 

1- E  Schematic  Drawing  of  Schlieren  System  60 

2- E  Schematic  Drawing  of  Coincidence 

Schlieren  System  64 

3- E  Schematic  Drawing  of  Coincidence 

Schlieren  System  65 


v 


Figure 

4-E 

1- F 

2- F 

3- F 

1-G 


Title  Page 

Schematic  Drawing  of  Schlieren  System 

with  Parallel  Light  Beams  Crossing 

Tunnel  Twice  67 

Isoclinic  Field  for  Determination  of 

Mach  Number  at  Exit  of  Segment  Two  of 

Nozzle  While  Operating  Shock  Free  89 

Isoclinic  Field  for  Determination  of 

Minimum  Diffuser  Height  89 

Isoclinic  Field  for  Determination  of 

Mach  Number  at  Minimum  Area  Section  89 

Decrease  of  Stagnation  Temperature  vs. 

Time  of  Tunnel  Operation  96 


vi 


LIST  OF  TABLES 


Table  Page 

I.  Design  Data  for  Mach  2.37  Nozzle  Blocks  40 

II.  Design  Data  for  Mach  1.50  Nozzle  Blocks  40 

III.  Data  for  Diffuser  set  at  Minimum  Starting 
Height  of  1.40  Inches  with  Mach  2.37 

Nozzle  Blocks  40 

IV.  Data  for  Diffuser  Set  at  Maximum  Height  of 

1.883  Inches  with  Mach  2.37  Nozzle  Blocks  40 

V.  Data  Sheet  for  Shock  Free  Operation  with 
Mach  2.37  Nozzle  Blocks,  Constant  Diffuser 

Area  and  Minimum  Back  Pressure  41 

VI.  Data  Sheet  for  Shock  Free  Operation  with 
Mach  2.37  Nozzle  Blocks,  Constant  Diffuser 

Area  and  Increased  Back  Pressure  41 

VII .  Data  Sheet  for  Shock  Free  Operation  with 

Mach  2.37  Nozzle  Blocks  and  Maximum  Dif¬ 
fuser  Area  41 

VIII.  Data  Sheet  for  Shock  Free  Operation  with 

Mach  2.37  Nozzle  Blocks  and  Minimum  Dif¬ 
fuser  Starting  Area.  41 

IX.  Data  Sheet  for  Shock-in-Nozzle  Operation 

with  Mach  2.37  Nozzle  Blocks  and  Constant 
Diffuser  Area  42 

X.  Simplified  Nozzle  Coordinates  Shock  Free 

Operation  90 

XI.  Simplified  Nozzle  Coordinates  for  Shock  in 

Nozzle  Operation  90 

XII.  Isoclinic  Data  for  Determination  of  Mach 
Number  in  Segment  Two  of  Nozzle  Operating 

Shock  Free  91 

XIII.  Diffuser  Contours  Required  to  Produce  Con¬ 
stant  Mach  Number  in  Test  Section  91 

XIV.  Isoclinic  Data  for  Determination  of  Min¬ 
imum  Diffuser  Area  92 

XV.  Isoclinic  Data  for  Determination  of  Mach 

Number  at  Minimum  Diffuser  Area  92 

XVI.  Summary  of  Adiabatic  Mach  Number  in  Nozzle  77 

vii 


LIST  OF  SYMBOLS 


cK 

Acute  Angle  between  Mach  Wavelet  and  horizontal 
line,  degrees 

A 

Area,  sq.  in. 

A* 

Area  in  tunnel  where  M  =  1.00,  sq.  in. 

*; 

Area  in  tunnel  upstream  of  shock  where 

M  =  1.00,  sq.  in. 

*y 

Area  in  tunnel  downstream  of  shock  where 

M  =  1.00,  sq.  in. 

CD 

Coefficient  of  discharge 

4) 

m 

Mass  rate  of  flow  -  lbm/sec. 

P 

Tunnel  pressure,  in.  Hg.a. 

A  P 

Differential  between  stagnation  pressure  and 
tunnel  pressure,  in.  Hg. 

^Porf 

Differential  pressure  across  orifice  plate, 
inches  of  water 

Pob 

Tunnel  back  pressure,  in.  Hg.a. 

Pon 

Stagnation  pressure  of  any  point  n  downstream 
of  shocking  system,  in.  Hg.  a. 

Pol 

Stilling  tank  stagnation  pressure,  in.  Hg.  a. 

P 

ox 

Stagnation  pressure  upstream  of  shock,  in.  Hg.a. 

p 

oy 

Stagnation  pressure  downstream  of  shock,  in.  hg.a 

T 

o 

Local  temperature,  R 

T  i 

ol 

Stilling  tank  stagnation  temperature,  °R 

viii 


INTRODUCTION 


There  are  two  basic  types  of  wind  tunnels  available  for  the  study 
of  compressible  fluid  flow  (1)  the  continuous  and  (2)  the  intermittent , 


With  continuous  wind  tunnels  the  problem  of  matching  com¬ 


pressor,  tunnel  volume  flow  and  pressure  ratio  characteristics  present 
major  difficulties  in  the  design  work.  In  addition  the  Mach  number  span 
of  a  continuous  wind  tunnel  is  usually  small  and  does  not  extend  from 
low  subsonic  velocities  to  high  supersonic  velocities.  These  factors 
combined  with  a  higher  cost  of  installation  and  operation  of  the  con¬ 
tinuous  tunnel  have  resulted  in  a  shift  towards  the  use  of  the  inter¬ 
mittent  tunnels. 

The  intermittent  tunnel  is  further  divided  into  two  basic  types^(l) 
the  vacuum  storage  drive  and  (2)  the  pressure  storage  drive,  Qll]]  •  In 
the  first  system,  atmospheric  air  passes  through  a  drier,  wind  tunnel 
and  diffuses  into  a  vacuum  vessel  which  is  continuously  evacuated  to  the 
atmosphere.  In  the  second  system,  dry  air  is  stored  under  pressure  in 
a  storage  tank;  the  air  passes  through  a  throttle  valve,  wind  tunnel 
and  exhausts  through  a  tunnel  diffuser  to  the  atmosphere. 

It  was  the  objective  of  this  thesis  to  install  and  evaluate  the  per¬ 
formance  of  a  small  pressure  storage  type  supersonic  wind  tunnel  system. 
The  components  available  at  the  start  of  the  project  consisted  of;  high 
pressure  compressor;  storage  tanks;  two  regulator  valves;  tunnel  test 
section;  and  Schlieren  system  components;  all  other  pieces  for  the 
system  had  to  be  either  specified  or  designed  and  specially  constructed. 
It  was  further  necessary  to  carry  out  assembly  and  installation  of  these 
components.  Having  completed  this  phase  of  the  objective,  it  was  then 

*  Numbers  in  brackets  refer  to  the  Bibliography  on  page  20. 
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possible  to  undertake  a  performance  evaluation  program.  This  con¬ 
sisted  of  both  a  qualitative  and  an  quantative  approach.  The  qualita¬ 
tive  evaluation  was  derived  primarily  from  Schlieren  photographs,  while 
the  quantitative  results  came  primarily  from  pressure  surveys. 

Figure  1  is  a  schematic  diagram  of  the  system  showing  the  basic 
design  and  necessary  components.  Discussion  of  the  design  and  per¬ 
tinent  details  of  certain  sections  of  the  system  is  found  in  Part  I. 

Figure  2  is  a  subassembly  drawing  of  the  small  tunnel  test  section 
for  which  the  performance  tests  were  made.  Results  of  these  tests  are 
presented  in  photographs,  graphical  form,  and  tabular  form.  Discussion 
of  these  tests  and  results  appear  in  Part  II. 

The  Appendicies  include  various  sample  calculations,  a  discussion 
of  Schlieren  systems,  operational  instructions,  design  of  a  differen¬ 
tial  pressure  transducer,  and  design  of  needed  heat  capacitors  for  the 
air  receivers. 
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Part  I.  Tunnel  System  Design 

A.  Air  Supply. 

Air  for  the  system  was  supplied  by  a  Chicago  Pneumatic  two  stage 
inner  cooled  air  compressor  rated  at  114  cubic  feet  of  free  air  per  min¬ 
ute  or  0.142  pounds  per  second.  The  compressor  is  driven  by  a  50  horse¬ 
power  induction  electric  motor.  The  compressor  discharge  pressure  can 

be  controlled  up  to  500  psi  by  the  regulator  valve  (Bailey  Type-XSDS-9M) . 

% 

Operational  instructions  for  the  compressor  are  given  in  Appendix  A. 

Air  from  the  compressor  goes  to  an  after-cooler  (Graham  Heliflow 
Type  9AA  145)  where  the  temperature  is  reduced  to  about  60°F,  and  then 
to  a  cyclone  type  separator  with  trap  (Wright  Austin),  where  the  con¬ 
densed  water  is  removed.  Calculations  show  the  moisure  content  after 
the  air  leaves  the  separator  to  be  about  0.000319  lb.  water  per  lb  of 
dry  air,  Appendix  B.  Drying  of  the  air  is  necessary  to  prevent  condensa¬ 
tion  from  occuring  in  the  tunnel  test  section,  ,  \\l~]  ,  . 

From  the  separator,  the  air  passes  through  an  oil  filter  (C.  M. 

Kemp,  Type  FPF-1610),  using  activated  alumina  (Vapoilsorb)  as  a  desi¬ 
ccant.  This  filter  was  modified  by  putting  a  60  mesh  stainless  steel 
strainer  on  the  inlet  side  along  with  a  fiberpacked  strainer.  This  modi¬ 
fication  was  made  to  prevent  blowback  of  desiccant  particles  to  the 
compressor  by  high  pressure  air  when  the  compressor  was  unloaded. 

The  air  then  goes  through  the  regulator  valve,  which  maintains 
a  constant  compressor  discharge  pressure,  to  the  two  air  receivers. 

The  receivers,  connected  in  parallel,  provide  a  total  storage  volume  of 
234  cubic  feet.  Figure  3  is  a  photograph  showing  these  various  compon¬ 
ents  in  relation  to  the  acoustic  enclosure  which  houses  the  compressor. 


3 


When  operating  the  wind  tunnel,  the  mass  flow  of  air  from  the 
receivers  to  the  tunnel  section  is  controlled  by  a  regulator  valve 
(Bailey,  Type  XSDS-9M).  The  regulator  valve  provides  for  stagnation 
pressure  control  up  to  200  psi.  Figure  4  shows  the  position  and  loca¬ 
tion  of  this  valve* 

B.  Tunnel  Section* 

The  tunnel  section  consists  of  six  basic  parts:  stilling  tank; 
entrance  nozzle;  test  section;  diffuser;  flow  meter;  and  back  pressure 
control* 

The  stilling  tank  was  made  from  a  30  inch  piece  of  nominal  6  inch 
pipe  flanged  at  both  ends.  The  flanges  were  faced  and  grooved  for  use 
with  0-ring  seals.  The  upstream  end  of  the  tank  was  connected  by 
means  of  a  piece  of  high  pressure  rubber  hose  to  the  regulator  valve 
discharge.  Figure  5  is  a  photograph  of  these  pieces.  In  the  upstream 
end  of  the  stilling  tank,  3  stainless  steel  60  mesh  screens  were  used 
to  diffuse  any  high  velocity  profile. 

The  entrance  nozzle  was  constructed  of  a  polyester  resin  and 
chopped  fiber  glass  molded  over  a  wooden  core.  A  pair  of  cubic  equa¬ 
tions  were  used  to  define  the  nozzle  contour  so  that  flow  separation 
would  not  occur.  The  flanges  were  molded  integral  with  the  nozzle  and 
were  faced  and  drilled  after  casting.  The  upstream  flange  has  a  1/2 
inch  lip  for  a  matching  connection  with  the  stilling  tank. 

The  test  section  was  designed  and  constructed  by  the  Ames 
Aeronautical  Laboratory,  Moffett  Field,  California.  Features  of  this 
test  section  consist  of  a  variable  area  supersonic  diffuser,  inter¬ 
changeable  nozzle  blocks,  and  interchangeable  side  walls  for  either 
visual  or  pressure  surveys. 
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The  present  sets  of  nozzle  blocks  are  for  Mach  numbers  1*50  and 
2.37,  The  sidewall  for  pressure  surveys  is  provided  with  32  pressure 
taps.  Figure  2  is  a  subassembly  drawing  of  the  section  and  Figure  6 
is  a  photograph  showing  the  section  with  side  walls  removed.  Tables  X 
and  II  give  the  nozzle  coordinates,  area  ratios,  and  design  Mach  numbers 
and  pressure  ratios  based  on  one-dimensional  calculations  for  the  two 
sets  of  nozzle  blocks. 

The  diffuser,  like  the  entrance  nozzle,  was  also  molded  over  a 
wooden  core  with  a  polyester  resin  and  fiber  glass.  To  provide  added 
strength,  1/8  by  3/8  inch  steel  bars  were  molded  into  the  side  and 
flanges.  The  design  was  straight  wall  with  a  maximum  divergance  angle 
of  7°  to  prevent  flow  separation;  and  diffusion  of  the  air  stream  was 
carried  to  a  point  where  the  flow  could  be  treated  as  incompressible  for 
metering  purposes. 

The  four  pieces  just  described  were  integrally  mounted  to  a  movable 
stand  and  are  shown  in  Figure  7. 

The  flow  meter  was  designed  according  to  A.S.M. E.  standards  for  a 
sharp  edge  orifice  with  D  and  (1/2)D  pressure  taps,  M  .  [?]  .  An  en¬ 
trance  length  of  20  pipe  diameters  and  an  exit  length  of  10  pipe  diameters 
were  provided.  Nominal  3  inch  pipe  was  used  which  was  bored  to  a  3.120 
inch  diameter  for  16  inches  either  side  of  the  orifice.  0-ring  seals 
were  used  throughout  including  the  flanged  joint  with  the  diffuser. 

Orifice  plates  were  made  of  1/8  inch  stainless  steel.  Appendix  C  presents 
the  details  for  computing  mass  flow  rates. 

The  back  pressure  control  consisted  of  a  manually  operated  3  inch 
globe  valve.  This  control  was  incorporated  so  that  changes  could  be  made 
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in  the  test  section  air  density.  Since  this  throttling  was  somewhat 
noisy,  a  large  truck  muffler  was  added  to  the  valve  exit. 

These  pieces  were  also  mounted  on  a  movable  stand  and  are  shown 
in  Figure  8. 

The  completely  assembled  tunnel  section  is  shown  in  the  photo¬ 
graph,  Figure  9. 

C.  Instrumentation. 

Instrumentation  for  the  tunnel  section  consisted  of  devices  for 
three  types  of  measurement;  pressure,  temperature,  and  velocity. 

Manometers  were  used  for  all  quantitative  pressure  measurements 
while  bourdon  tube  gages  were  hsed  for  qualitative  check  purposes. 

Two  120  inch  mercury  differential  manometers  were  utilized  with  the 
32  tap  pressure  survey  sidewall  of  the  test  section.  Two  16  connector 
manifolds  with  quick  opening  valves  were  made  to  accommodate  the  32  taps 
to  the  two  manometers.  Figure  10  shows  the  arrangement  of  the  manifolds 
with  the  manometers.  It  was  originally  planned  to  produce  strain  gage 
type  differential  pressure  transducers  for  these  measurements  to  de¬ 
crease  the  reading  time  required;  fabrication  difficulties,  however, 
prevented  the  completion  of  this  project.  Appendix  D. 

A  120  inch  mercury  differential  manometer  was  also  used  to  monitor 
the  test  section  throat  pressure  continuously  during  pressure  survey  runs. 

Stagnation  pressure,  from  the  stilling  tank  tap,  was  likewise  con¬ 
tinuously  monitored  with  a  120  inch  mercury  manometer  besides  being  fed 
to  the  differential  manometers. 

Static  pressure  for  the  flow  meter  could  be  read  with  either  an 
oil  (sp.  gr,  =  0.834)  or  mercury  filled  60  Inch  manometer.  Differen¬ 
tial  pressure  for  the  orifice  was  determined  with  a  water  manometer. 
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Ellison  inclined  manometers  were  used  in  conjunction  with  the  pitot 

tube. 

The  only  quantitative  temperature  measured  was  the  stagnation 
temperature  in  the  stilling  tank.  This  was  accomplished  with  a 
stainless  steel  radiation  shielded  30  gage  copper-constantan  thermocouple 
referenced  to  ice.  The  output  emf  of  the  thermocouple  was  read  with  a 
continuous  balance  precision  potentiometer  (Minneapolis-Honeywell8  Model 
SY156X-VH).  Various  other  temperatures  of  the  tunnel  system  were  meas¬ 
ured  with  gas  bulb  indicating  thermometers. 

Velocity  measurements  in  the  stilling  tank  were  made  with  a  1/8 
inch  diameter  Prandt  1-Pitot  probe.  These  measurements  were  made  sole¬ 
ly  for  the  purpose  of  determining  the  uniformity  of  the  velocity  pro¬ 
file  leaving  the  straightening  screens  in  the  stilling  tank  and  were 
only  of  a  qualitative  nature.  The  Pitot  probe  was  removed  from  the 
stilling  tank  after  finding  that  the  screens  adequately  diffused  the 
high  velocity  entering  air  to  a  uniform  low  velocity  profile. 

Velocity  measurements  and  visual  flow  studies  in  the  test  sec¬ 
tion  were  accomplished  with  a  coincidence  type  Schlieren  systeii  manu¬ 
factured  by  Tinsley  Laboratories  Inc.  A  schematic  diagram  is  shown  in 
Figure  11.  The  components  were  mounted  to  a  rigid  moveable  stand  with 
the  wheels  set  in  tracks  parallel  to  the  test  section.  Figure  8<>  The 
viewing  section  of  the  Schlieren  system  was  modified  by  installation 
of  an  Ilex  Synchro  shutter  and  a  4  by  5  inch  ground  glass  view  plate 
and  film  pack  holder.  Figure  12.  This  modification  allowed  photographs 
of  the  flow  behavior  to  be  made  which,  when  enlarged  in  print ing^ pro¬ 
vided  for  qualitative  and  some  quantitative  measurements.  Appendix  E 


7 


presents  a  general  theory  of  the  Schlieren  method  of  flow  visualization 


as  well  as 


operating  instructions  for  this  particular  system. 
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Part  II.  Performance  Evaluation 


A.  Qualitative  Results. 

The  qualitative  evaluation  of  the  tunnel  system  performance  was 
derived  primarily  from  Schlieren  photographs.  Figure  13  depicts  the 
tunnel  section  operation  under  shock  free  conditions.  The  knife  edge 
intercepted  the  light  beam  in  both  a  vertical  and  horizontal  position 
to  show  the  change  of  density  patterns  in  the  horizontal  and  vertical 
directions.  Figure  13c,  with  the  knife  edge  horizontal,  clearly  shows 
the  boundary  layer. 

Criss-crossing  the  tunnel  outline  are  Mach  wavelets  produced  by 
wall  friction  and  other  irregularities  of  the  tunnel  surfaces.  These 


wavelets  are  related 


to  the  Mach  number  of  the  air  flow  in  the 


tunnel  by  the  relationship, 

•  ^  -  1 
sin  0<  -  pj 

where  <?C  is  the  acute  angle  between  the  wavelets  and  a  horizontal  line. 
Figure  18a  is  an  enlargement  of  a  section  of  Figure  13a  and  Figure  18b 
is  an  enlargement  of  Figure  13b.  The  angle  c/\  in  Figure  18a  is  27.50° 
with  a  resulting  Mach  number  of  2,17.  This  wavelet  originated  at  the 
junction  of  the  nozzle  and  the  test  section.  Utilizing  the  pressure 
ratios,  the  Mach  number  at  the  corresponding  point  was  determined  to 
be  2.125.  The  angle  in  Figure  18b  is  29.25°  with  a  resulting  Mach 
number  of  2,04.  The  Mach  number  calculated  for  this  position  was  1.828, 
This  wavelet  originated  at  approximately  pressure  tap  30. 

Figure  14  depicts  a  normal  shock  occuring  in  the  diffuser  section 
which  was  being  held  at  a  constant  height  of  1.718  inches.  Inner  action 
between  the  shock  and  the  boundary  layer  is  clearly  indicated.  Up¬ 
stream  of  the  shock  the  boundary  layer  is  stable.  At  station  12  there 
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is  a  separation  of  the  boundary  layer  due  to  the  interaction  of  the 
shock  and  the  boundary  layer.  Less  severe  interaction  occurs  between 
the  boundary  layer  and  the  repeated  shocks  downstream  until  complete 
subsonic  flow  is  established  downstream  of  station  15. 

Comparing  Figure  14a  operating  under  a  stagnation  pressure  of 
67  inches  of  mercury  and  a  back  pressure  of  30.5  inches  of  mercury  to 
Figure  14c  for  tunnel  operation  with  a  stagnation  pressure  of  81.5  inches 
of  mercury  and  a  back  pressure  of  38.3  inches  of  mercury,  it  may  be  seen 
that  the  resultant  shock  has  moved  upstream  almost  three  inches.  The 
pressure  ratio  at  station  11  in  Figure  14a  where  no  shock  exists  is 
approximately  0.1088  as  taken  from  Figure  23  or  Table  V.  At  a  point 
forward  of  the  shock  in  Figure  14c  a  pressure  ratio  of  0.103  exists 
as  taken  from  Table  VII. The  increased  density  of  the  air  and  the  high¬ 
er  Reynolds  number  must  account  for  the  shock  moving  upstream. 

Reduction  of  the  stagnation  pressure  also  moves  the  shock  up¬ 
stream  as  shown  in  the  sequence  of  photographs  in  Figure  15.  A  strong 
shock  exists  at  the  nozzle  exit  of  Figure  15a  with  a  stagnation  pressure 
of  65.35  inches  of  mercury.  Reduction  of  the  stagnation  pressure  to 
59.96  inches  of  mercury  allowed  the  shock  to  move  upstream  one  and  one 
half  inches  towards  the  nozzle  throat  as  shown  in  Figure  15c.  Reduction 
of  the  shock  intensity  may  be  noted  also  as  the  shock  moves  upstream  by 
comparing  Figures  15a  and  15c. 

The  minimum  starting  area  of  the  variable  diffuser  was  found  to 
range  between  0.695  and  0.700  square  inches  as  compared  to  the  isentropic 
minimum  area  of  0.665  square  inches.  The  area  actually  required  corres¬ 
ponds  to  a  diffuser  heights  of  1.39  to  1.40  inches.  Figure  16a  is  a 
photograph  of  the  diffuser  nozzle  set  at  a  height  of  1.39  inches.  The 
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shock  could  not  be  blown  through  the  diffuser  with  a  stagnation  pressure 

of  82.71  inches  of  mercury.  Increasing  the  diffuser  height  to  1,40 

inches,  (area  change  of  0.005  square  inches)  permitted  the  shock  to  be 

blown  through  the  diffuser  nozzle  and  the  stagnation  pressure  reduced 

* 

as  shown  in  Figure  16b.  The  stagnation  pressure  could  have  been  reduced 
to  move  the  shock  farther  upstream  to  the  throat  position.  Figure  17a 
and  17b  depict  tunnel  operation  with  the  diffuser  height  set  slightly 
greater  than  the  required  minimum  and  show  the  shock  standing  in  the 
diffuser  throat  and  the  stagnation  pressure  reduced  to  a  minimum  while 
maintaining  the  shock  location.  Several  minutes  after  Figure  17a  was 
taken  the  stagnation  pressure  was  reduced  from  60.00  inches  of  mercury  to 
59.50  inches  of  mercury  and  the  shock  moved  through  the  diffuser  nozzle 
and  relocated  itself  in  the  diverging  portion  of  the  tunnel  nozzle. 

The  sharpness  of  the  Schlieren  photographs  was  poor  due  to  the  use 
of  non-collimated  light  and  the  fact  the  light  beam  passed  through  the 
test  region  twice,  Appendix  E.  Figure  16c  gives  an  indication  of  the 
seriousness  of  this  problem.  The  knife  edge  was  not  inserted  into  the 
light  beam.  The  black  parallel  vertical  lines  are  images  produced  by 
placing  fine  copper  wires  on  the  outside  of  the  Schlieren  window  nearest 
to  the  spherical  mirror;  the  diagonal  parallel  lines  are  images  produced 
by  placing  the  same  size  wire  on  the  outside  of  the  Schlieren  window 
nearest  to  the  light  source.  The  greater  the  distance  between  the 
images  of  any  wire,  the  farther  the  light  beam  had  to  travel  before 
reaching  the  spherical  mirror.  The  non-parallel  light  beam  intercepted 
the  diagonal  wire  first  and  hence  the  image  of  this  wire  was  deflected 
more  due  to  the  longer  distance  to  the  spherical  mirror  than  it  was  for 
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any  of  the  vertical  wires  which  were  closer  to  the  spherical  mirror*  As 
evidenced  by  the  two  images  formed  for  each  wire*  the  Schlieren  picture 
of  the  of  the  air  flow  would  tend  to  be  fuzzy  due  to  superimposing  the 
two  images  formed.  The  horizontal  black  line  is  a  wire  placed  at  the 
tunnel  centerline. 

Also  contributing  to  poor  quality  of  the  Schlieren  pictures  was 

the  incomplete  removal  of  oil  vapor  from  the  air.  During  tunnel  opera¬ 
's 

tion,  oil  would  be  deposited  in  a  very  thin  film  on  the  glass  wall  and 
cause  a  refraction  of  the  light.  In  addition,  the  film  deposited  was 
not  of  a  uniform  nature.  This  condition  may  be  seen  in  Figure  17c. 

The  dark  convex  lines  at  stations  4  and  5  are  the  result  of  an  oil  film 
deposit.  These  same  lines  could  still  be  seen  after  securing  the  air 
supply  and  leaving  the  Schlieren  unit  in  operation.  Due  to  this  oil 
film  formation,  it  was  necessary  to  remove  the  Schlieren  window  at 
frequent  intervals  and  wipe  off  the  small  amount  of  oil  that  had  been 
deposited.  Water  vapor  in  the  atmosphere  condensing  on  the  outside 
surface  of  the  Schlieren  window  was  also  a  factor  in  obtaining  clear 
Schlieren  photographs. 

As  the  vertical  wires  were  spaced  at  one  inch  intervals,  the  dis¬ 
tance  between  wires  in  the  photograph  may  be  measured  and  comparison  of 
this  measured  value  to  the  known  one  inch  dimension  permits  determining 
a  scale  factor  for  the  photograph.  This  method  was  utilized  to  deter¬ 
mine  the  diffuser  heights  for  heights  other  than  1.718  inches.  A 
lucite  gage  block,  1.718  inches  high,  was  utilized  when  runs  of  this 
diffuser  height  were  made. 
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The  apparent  difference  in  intensity  of  the  shock  and 
the  visibility  of  the  boundary  layer  caused  by  the  knife 
edge  orientation  demonstrates  why  a  horizontal  knife  is 
generally  used  when  studying  boundary  layers  and  a  vertical 
knife  is  used  when  studying  shock  waves. 

B.  Quantitative  Results . 

The  quantitative  evaluation  was  derived  primarily  from 
pressure  surveys.  The  Mach  numbers  actually  obtained  in 
the  tunnel  were  generally  less  than  the  predicted  adiabatic 
Mach  numbers  due  to  friction  effects  not  accounted  for,  ob¬ 
lique  shocks  and  other  departures  from  a  one  dimensional 
flow  analysis.  Figure  20  is  a  plot  of  the  adiabatic  Mach 
numbers  obtained  from  the  pressure  and  area  ratios  for  a 
constant  diffuser  height  of  1.718  inches  vs.  distance  down¬ 
stream  of  the  tunnel  nozzle  throat.  The  solid  curve  is  the 
isentropic  Mach  number  based  on  the  change  of  area  ratios 
only  as  given  in  Table  I,  while  the  dashed  curve  is  the 
predicted  adiabatic  curve  as  determined  by  the  use  of  iso¬ 
clines  in  Appendix  F.  The  circles  are  the  adiabatic  Mach 
numbers  obtained  from  the  data  taken  on  the  tunnel  and  are 
tabulated  in  Table  V,  based  on  the  procedure  outlined  in 
Appendix  F. 

Up  to  a  Mach  number  of  1.936  at  pressure  tap  number 
4,  the  actual  Mach  numbers  followed  the  adiabatic  curve 
closely,  exceeding  both  the  isentropic  and  adiabatic 
curves  between  pressure  taps  numbers  4  and  7.  At  this 
point,  the  Mach  numbers  obtained  in  the  tunnel  began  to 
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fall  below  the  isentropic  curve  but  still  above  the  adia¬ 
batic  curve  to  pressure  tap  number  9.  The  maximum  Mach 
number  of  2.275  was  obtained  at  pressure  tap  number  7.  The 
decrease  in  Mach  number  at  taps  number  10  and  24  is  un¬ 
explained.  This  decrease  in  Mach  number  at  these  pressure 
taps  appeared  consistently  in  all  runs  made  during  the 
evaluation  of  the  tunnel. 

An  attempt  was  made  to  maintain  a  constant  Mach  number 
in  the  test  section  by  increasing  the  height  of  the  varia¬ 
ble  area  diffuser.  The  results  of  this  run  are  shown 
graphically  in  Figure  21  and  in  tabular  form  in  Table  VII. 
Because  of  design  limitations,  the  diffuser  height  could 
not  exceed  1.883  inches.  The  calculated  maximum  height  of 
the  diffuser  to  maintain  a  constant  Mach  number  in  the  test 
section  including  the  effects  of  friction  was  2.145  inches, 
Appendix  F,  based  on  a  friction  factor  of  4f  =  0.01.  Table 
IV  presents  the  diffuser  coordinates  vs.  diffuser  length  in 
tabular  form  while  Figure  19b  presents  the  same  data  in 
graphical  form.  Examination  of  Figure  21  indicates  a 
definite  drop  in  Mach  number  between  the  nozzle  exit  and 
tap  24,  in  much  the  same  manner  as  the  drop  for  constant 
diffuser  area.  Figure  20,  between  the  same  pressure  tap 
locations.  As  the  diffuser  did  not  increase  in  height 
appreciably  in  the  first  1.5  inches  of  diffuser  length,  the 
effect  of  friction  was  to  continually  reduce  the  Mach  num¬ 
ber.  Recovery  of  the  Mach  number  was  accomplished  shortly 
after  the  diffuser  height  began  to  increase.  The  recovered 
Mach  number  is  somewhat  higher  than  that  for  constant 
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diffuser  area,  and  appears  to  be  maintaining  a  more  con¬ 
stant  value,  thus  indicating  the  effect  of  increasing  the 
area  of  the  diffuser.  Had  it  been  possible  to  increase  the 
diffuser  height  to  the  calculated  height  of  2.145  inches, 
the  resulting  Mach  number  would  have  been  increased  much 
more.  Until  the  behavior  in  Mach  number  without  area 
change  can  be  explained,  the  effect  of  changing  the  area  of 
the  diffuser  can  not  be  adequately  predicted. 

Figure  24  is  a  plot  of  adiabatic  Mach  numbers  vs. 
distance  from  the  tunnel  nozzle  throat  for  two  operating 
conditions:  (1)  the  diffuser  height  set  at  constant  1.718 

inches  and  a  back  pressure  of  38.98  inches  of  mercury,  and 
(2)  the  diffuser  set  at  the  minimum  starting  height  of  1.40 
inches  with  a  back  pressure  of  30.54  inches  of  mercury. 

The  minimum  diffuser  height  was  based  on  the  actual  minimum 
height  required  to  operate  the  tunnel  supersonically  in  the 
test  section  confirmed  by  the  calculated  height  in  Appendix 
F.  The  calculated  minimum  height  for  isentropic  flow  was 
found  to  be  1.33  inches;  while  that  for  adiabatic  flow  was 
1.40  inches.  Appendix  F.  Table  III  gives  the  diffuser  co¬ 
ordinates  in  tabular  form  and  Figure  19b  presents  the  same 
information  graphically.  It  should  be  noted  there  was 
very  little  change  of  area  of  the  diffuser  near  the  nozzle 
outlet  when  adjustment  of  the  diffuser  height  was  made 
downstream. 

The  data  followed  the  predicted  curves  in  much  the 
same  general  pattern  as  when  the  tunnel  was  operating  under 
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a  constant  diffuser  height  shock  free  condition.  The 
characteristic  drop  in  Mach  numbers  may  be  noted  at  pres¬ 
sure  taps  10  and  24.  A  predicted  adiabatic  Mach  number  at 
the  diffuser  throat  was  calculated  with  the  diffuser  height 
set  at  1.40  inches.  Appendix  F.  This  predicted  value  of 
M  =  1.64  does  not  correlate  too  well  with  the  value  of 
M  =  1.34  obtained  in  the  tunnel  as  shown  in  Figure  24. 

While  operating  supersonically,  the  tunnel  had  a  mass 
flow  rate  ranging  from  0.25  lbm/sec.  to  0.36  lbm/sec. ,  as 
determined  by  the  use  of  an  orifice  meter.  The  correspond¬ 
ing  stagnation  pressure  varied  from  62  to  88  inches  of 
mercury.  The  mass  rate  of  flow  for  the  orifice  meter  may  be 
determined  by  the  procedure  outlined  in  Appendix  C,  or  the 
ideal  mass  rate  of  flow  for  the  tunnel  nozzle  throat,  as 
determined  by  Fliegner's  formula,  may  be  corrected  by  the 
use  of  the  empirical  Coefficient  of  Discharge  presented 
graphically  in  Figure  25. 

The  average  friction  factor  4f,  in  the  diffuser  sec¬ 
tion  was  found  to  be  0.01865  when  the  tunnel  was  operating 
shock  free  supersonically.  Appendix  F.  The  diffuser  height 
was  maintained  at  a  constant  1.718  inches  and  the  data 
from  Table  V  utilized.  All  predictions  made  involving  the 
use  of  friction  were  based  on  the  assumption  the  mean 
friction  coefficient  4f  was  0.01  and  was  constant  at  this 
value. 

Figures  22  and  23,  and  Table  IX  present  the  results 
of  a  pressure  traverse  of  the  tunnel  operating  with  a 


16 


normal  shock  in  the  tunnel  nozzle.  Figure  22  is  a  plot  of 
the  Mach  number  vs.  distance  from  the  nozzle  throat.  The 
solid  line  represents  the  isentropic  case  assuming  a  single 
normal  shock  occurring  at  M  =  2.05  and  the  dashed  line  the 
adiabatic  case  assuming  a  normal  shock  at  M  s  2.00.  Imme¬ 
diately  downstream  of  the  single  shock  the  isentropic  Mach 
number  will  be  decreased  at  0.569  and  continue  to  decrease 
in  the  nozzle,  due  to  the  increasing  area,  until  it  reaches 
a  value  of  0.388  at  the  end  of  the  nozzle;  the  isentropic 
Mach  number  then  remains  constant  in  the  diffuser  section. 
For  the  adiabatic  case,  the  Mach  number  immediately  down¬ 
stream  of  the  shock  will  be  0.577,  based  on  the  shock  occur¬ 
ring  at  M  »  2.00.  The  Mach  number  continues  to  decrease  for 
the  subsonic  flow  but  at  a  slower  rate.  Between  pressure 
taps  9  and  15,  the  area  change  is  slight  and  the  Mach  number 
remains  constant  at  M  -  0.429.  In  the  constant  area  dif¬ 
fuser,  the  adiabatic  Mach  number  is  increased  due  to  the 
frictional  effects.  As  shown  in  the  Schlieren  photographs 
by  the  repeated  shock  pattern.  Figures  13  to  18,  the  flow 
conditions  through  a  normal  shock  are  definitely  not  a 
simple  one  dimensional  process.  A  shock  occurs  and  the  flow 
immediately  goes  subsonic.  Due  to  inner  action  with  the 
boundary  layer,  the  flow  accelerates  to  a  lower  supersonic 
value  and  the  process  is  repeated  until  the  flow  remains 
subsonic.  The  overall  effect  of  Mach  number  reduction 
caused  by  the  multiple  shocks  is  shown  in  Figure  22.  Due 
to  the  complex  inner  action  of  the  normal  shocks  and  the 
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oblique  shocks  with  the  boundary  layer,  and  the  change  of 
stagnation  pressure  for  each  shock  occurring,  definite 
determination  of  the  Mach  number  through  a  shocking  system 
is  difficult.  The  method  for  calculating  the  Mach  numbers 
of  Figure  22  is  presented  in  Appendix  F. 

At  a  Mach  number  exceeding  1.00  the  effect  of  friction 
is  to  reduce  the  Mach  number  to  a  limiting  value  of  1.00. 

At  Mach  numbers  lower  than  unity,  the  effect  of  friction  is 
reversed;  namely,  the  subsonic  flow  is  accelerated  to  a 
limiting  Mach  number  of  loOO.  This  increase  of  Mach  number 
due  to  the  effects  of  friction  on  subsonic  flow  may  be  noted 
by  the  gradual  increase  of  Mach  number  of  the  subsonic  flow 
in  Figure  22  or  in  Table  IX. 

Figure  23  is  a  plot  of  the  tunnel  pressure  ratios 
based  on  the  stilling  tank  stagnation  pressure  vs.  distance 
from  the  nozzle  throat  for  the  cases  of  shock  free  and 
shock  in  the  nozzle  operations.  The  pressure  ratio  pattern 
for  shock  free  operation  follows  the  isentropic  curve 
closely  with  deviations  commencing  about  four  inches  from 
the  nozzle  throat.  With  shock  in  the  nozzle  operation, 
failure  of  the  pressure  ratio  downstream  of  the  shock  to 
recover  to  the  isentropic  value  is  a  loss  of  stagnation 
pressure  caused  by  the  repeated  shock  pattern. 
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Conclusions  and  Recommendations 

A.  Conclusions. 

1.  The  qualitative  and  quantitative  evaluation  of  the  tun¬ 
nel  system  as  installed,  show  it  to  be  a  suitable  laboratory 
device  for  the  study  of  compressible  fluid  flow. 

2.  The  minimum  starting  height  of  the  variable  area  dif¬ 
fuser  was  found  to  be  1.40  inches. 

3.  The  average  friction  factor,  4f,  was  found  to  be  0.02044, 
but  is  believed  high  in  value. 

4.  Schlieren  photographs  of  the  flow  field  may  be  taken  in 
black  and  white  or  in  color. 

5.  A  vertical  knife  edge  orientation  in  the  Schlieren  sys¬ 
tem  is  best  for  studying  shock  phenomenon. 

6.  A  horizontal  knife  edge  orientation  in  the  Schlieren 
system  is  best  for  studying  the  boundary  layer. 

7.  The  sharpness  of  the  Schlieren  photographs  is  poor  due 
to  the  use  of  non-collimated  light  and  the  light  beam  pass¬ 
ing  through  the  tunnel  test  section  twice. 

8.  The  phenomenon  of  a  standing  shock  is  definitely  not 
one  dimensional. 

B.  Recommendations . 

It  is  recommended  the  following  modifications  for  im¬ 
proving  the  performance  and  usefulness  of  the  tunnel  system 
be  undertaken : 

1.  An  additional  relay  be  installed  on  the  Bailey  stagna¬ 
tion  pressure  controller  to  maintain  a  more  uniform  stag¬ 
nation  pressure. 
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2.  Additional  heat  capacitors  be  installed  in  the  receiv¬ 
ers  to  maintain  a  more  constant  stagnation  temperature. 

3.  A  total  pressure  probe  be  installed  in  the  fixed  geom¬ 
etry  diffuser  to  obtain  the  back  pressure  stagnation  pressure. 

4.  For  more  quantitative  visual  results,  the  Schlieren 
system  be  modified  to  a  single  pass  parallel  light  beam 
system. 

5.  A  more  rapid  means  of  obtaining  the  differential  pres¬ 
sure  measurements  be  installed. 

6.  The  camera  system  be  modified  to  a  reflex  type  system. 

7.  An  additional  oil  filter  be  installed. 

8.  Additional  pressure  taps  be  installed  downstream  of  the 
variable  area  diffuser  in  the  existing  pressure  survey  tun¬ 
nel  sidewall. 

9.  The  geometry  of  the  variable  area  diffuser  be  modified 
to  permit  constant  Mach  number  runs. 
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Schematic  Diagram  of  the  Tunnel  System 
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Figure  5.  Stilling  Tank. 


Figure  10.  Manometer  Manifolds 
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Figure  6.  Test  Section,  Disassembled 
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Figure  7.  Stilling  Tank,  Entrance  Nozzle,  Test  Section 
and  Diffuser  Assembly. 
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Figure  9.  Overall  View  of  Tunnel  Section  from  Front 
Side. 


Figure  8.  Overall  View  of  Tunnel  Section  from  Rear 
Side. 
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FIG.  II  .  3 chueren  Camera  System 


Figure  12.  Schlieren  System  Light  Source,  Knife  Edge, 
and  Camera -Viewing  Components. 


A.  Vertical  Knife.  P0  -  71.46  In  Hg.  a.,  Ph  -  30.92  to.  Hg.  a.,  TQ  -  61.1*  P. 
Korda  Running  Shock  rraa  with  Air  Flnv  Right  to  Left.  Kota  Strong  Mach 
Uavalat  Originating  at  Junction  of  Nozzle  and  Taat  Sactlnn.  Saa  Flgura  18. 


B.  Vertical  Knife.  Pn  ■  72.31  In  Hg.  a.,  Pj,  »  30.92  In.  Hg.  a.  f  T0  •  59.23  F. 
Tunnel  Running  Shock  Free  With  Air  Flnv  Right  to  Left. 


C.  Horizontal  Knife.  P0  -  72.71  to.  Hg.  a.,  Pv  -  30.92  in.  Hg.  a.,  Tn  •  65.5*  F. 
Tunnel  Running  Shock  Fraa  With  Air  Flnv  Right  tn  Left.  Kota  Thtcknasa  of  tha 
Boundary  Layers. 


Tunnel  Running  Shock  Fraa 


Jc 


Flgura  13. 


A.  Vertical  Knife, 
Dlffuaer  Haight 


P  -  67.35  In.  Hg.  a.,  Pb  -  30.45  In.  Hg.  a.,  T  =  66.3*  P. 
1.718  lnchea.  Normal  Shock  at  Station  12. 


B.  Horizontal  Emit* 
Dlffuaer  Haight 


F  •  67.41  In.  Hg.  a.,  Pb  -  30.86  In.  Hg. 
1.718  Inchaa.  Normal  Shock  at  Station  12. 


T0  -  65.9*  F. 


C.  Vertical  Knife.  PQ  *  81.52  In.  Hg.  a.,  Pfe  -  38.32  In.  Hg.  a.,  T  -  57.0*  P. 
Dlffuaer  Height  1.718  Inchaa.  Due  to  Higher  Danalty  and  Raynolde  Number 
Shock  haa  Moved  Upatream  aa  Coopered  tn  Figure  14A. 


Figure  14.  Effect  of  Incraaaed  Back  Praaaure  nn  Shock  Location 
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A.  Vertical  Knife.  PQ  .  65.35  In.  Hg.  a.,  Pb  .  30.65  In.  Hg.  a.,  Tn  -  59.6*  P. 
Boundary  Layer  Vlalbla  on  Nozzle  Walla  and  Incraaatng  In  Thlckneee  Dovnatraen 
of  Shock. 


5lfV  Po  “  ln*  •*.  ?b  “  30*83  In.  Hg.  a.,  T0  -  63.2*  P. 

Note  Pull  off  of  Boundary  Layer  at  Shock  and  Aeymetrlcal  Shading  of  Cnrraa- 
pnndlng  parte  on  Oppnelte  Sldaa  of  Tunnel  Centerline.  Shock  la  Moving  Tnvarde 
Nozzle  Throat  aa  Prt  la  Decraaaad  fro*  Figure  15  A. 


C.  Vertical  Knlfa.  P„  -  59.96  lo. 
Weaker  Shock  Moving  Opetreen  aa 


Hg.  a.,  Pb  -  30.73  In.  Hg.  a.,  Tn  -  61.8* 
Pn  la  Reduced  an  Additional  Aatount. 


P. 


Figure  15.  Effect  nf  Reducing  Stagnatloo  Praaaura  for  Shock  lo  Nnzzla 
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A.  Vertical  Knife. 
Shock  Could  not 


PQ  -  82.71  In.  Hg . a , ,  T0  -  57.6* 
ba  Blown  Through  Diffuser  Nozzle. 


F,  Dlffuaar  Haight 


B.  Vertical  Knife.  Pc  -  60.26  in.  Hg.  a.,  T0  -  58.6*  F.  Dlffoaar  Haight  -  1.40  lo. 
Area  hea  bean  Incraaaad  0.005  Square  tochee  Compared  to  Figure  16  A,  aod  Shock 
haa  baen  Blown  Through  Dlffuaar  Nozzle.  Praaaora  Coold  ba  Reduced  more  to  we 
Shock  towards  the  Throat. 


C.  Tunnel  oot  Running.  Knife  Edge  not  Iosartad.  Effect  of  Noo-Parellel  Two  Peas 
Light  la  Evident.  Double  Vartlcal  Lloaa  are  Image*  of  Wlree  Placed  oo 
Spherical  Mirror  Side  oS. Schllaran  Window.  Diagonal  Lloaa  are  Images  of  Wlraa 
Placed  on  Light  Source  Side  of  Schllaran  Window.  Light  Paselng  These  Wlraa 
had  to  Travel  Farther  before  .Striking  Spherical  Mirror  and  Hence  Mora  Deflection 
of  the  Imaee. 


Figure  16.  Effect  of  Area  Change  on  Shock  Location  In  Dlffuaar  Noztie. 
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A.  Vartlcal  Kolfa.  P0  -  60.00  in.  Rg.  T0  -  55.0*  P.  Dlffuaar 
Haight  •  1.42  inch** .  Shock  Pa  Had  Back  Through  Difftutr  goes  la 
ahao  PA  mi  Dacraaiad  to  59. SO  In.  Rg.  a. 


>.  Horizontal  tnifa.  P0  •  60.61  io.  Rg.  a.,  T0  •  63.1*  F.  Dlffuaar 
Haight  -  1.47  lachaa.  Rota  Strong  Pull  Off  of  Bowidarjr  Layar  at 
Station  12.25. 


C.  Vart teal  Knifa.  P  •  72.31  In.  Hg.  a.,  T0  •  5S.7*  F.  Dark  Coocara 
Lina  a  at  Statiooa  4  and  5  an  tha  Raault  of  011  Fllu  Dapoalt  oo  Walla. 
Dark  Barltontal  Lina a  Starting  at  Station  3  urn  alao  oil  Straaka  on 
Schllaraa  Walla. 


Flgura  17.  Tunaal  ftaaalag  With  Radacad  Dlffoaar  Araa. 
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A.  Vertical  Knife.  PQ  -  71.46  lo.  Hg.  e.  Pb  -  30.92  lo.  Hg.  e  T  -  61.1*  F. 
Air  Flow  U  tight  to  Left.  Section  le  eo  Eolargeaant  of  Figure  13  A. 

•*®**1*  1*  Fanning  Shock  Free.  Mach  Rtnber  datermload  fro®  Wevelete  by 
Ralat looeblp  S m  <K  *  .At  Corraepoodlng  Polot  Adlebetlc  Mach 

Ntaibtr  le  2.125.  ^ 


B.  Vertical  Knife.  PQ  -  72.31  lo.  Hg.  a.,  Pfc  -  30.92  lo.  Hg.  a.,  Te  -  59.23*  F. 
Air  Flow  la  Right  to  Left.  Section  le  an  Eolargemaot  of  Flgare  13B.  At 
Corraepoodlng  Point  Adiabatic  Mach  (twber  la  1.828. 


Figure  18.  Determlnat loo  of  Mach  Number  from  Mach  Wevelete. 
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A.  Vertical  Knlfa.  PQ  -  71.46  to.  Hg.  a.  Pfe  -  30.92  to.  Hg.  a.  T  -  61.1*  F. 
Air  Flo*  ta  Fight  to  Left.  Section  la  eo  Eolergament  of  Figure  13  A. 

Nora  la  la  Running  Shock  Free.  Mach  Number  determined  from  Wavelets  by 
Relationship  S //?  e*  =  .At  Corresponding  Point  Adiabatic  Mach 

Ncaber  la  2.125. 


B.  Vertical  Inlfa.  PQ  -  72.31  In.  Hg.  a.,  Pb  -  30.92  In.  Hg.  a.,  T0  -  59.23*  F. 
Air  Flow  Is  Right  to  Left.  Sectloo  la  an  Enlargement  of  Figure  13B.  At 
Corraapondlng  Point  Adiabatic  Mach  Niabsr  la  1.828. 


Figure  18.  Datarmlnatlon  of  Mach  Number  from  Mach  Wevalete. 
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APPENDIX  A 


TUNNEL  SYSTEM  OPERATION  PROCEDURE 

I.  For  charging  the  storage  tanks  the  following  procedure  should  be 
followed: 

1.  Open  the  3/4  inch  valve  connecting  the  trap  with  the  separa¬ 
tor. 

2.  Close  the  3/4  inch  drain  valve  located  below  the  tee  after 
the  trap. 

3.  Open  the  1/4  inch  balance  line  valve  between  the  separator 
and  the  trap. 

4.  Open  the  cooling  water  valve  to  the  after-cooler  heat  ex¬ 
changer. 

5.  Close  both  two  inch  inlet  valves  located  at  the  entrance  to 
the  low  pressure  storage  tanks. 

6.  Open  the  valve  on  the  low  pressure  after-cooler. 

NOTE:  Items  5  and  6  are  safety  precautions  to  prevent  high 
pressure  air  from  entering  the  low  pressure  tanks. 

7.  Close  the  drain  valve  on  each  high  pressure  tank. 

8.  Close  the  two  inch  discharge  valve  on  each  high  pressure 
tank  outlet. 

9.  Open  the  inlet  valve  to  each  high  pressure  tank. 

10.  Place  the  unloader  on  the  compressor  in  the  START  position 
(vertical  with  the  index  head  down.) 

11.  Place  cover  on  the  compressor  oil  vent. 

12.  Open  the  air  supply  valve  to  the  Bailey  controller  (Pressure 
should  read  28  psi). 
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13.  Adjust  compressor  discharge  pressure  for  a  pressure  slight = 
ly  higher  than  any  pressure  indicated  in  the  receivers.  This 
is  done  by  turning  the  adjustment  knob  on  the  Bailey  Con¬ 
troller  counterclockwise  until  the  red  hand  indicates  the  de¬ 
sired  pressure. 

14.  Connect  the  cooling  air  circulating  centrifugal  blower  in 
the  compressor  shack. 

15.  Throw  the  master  power  switch  located  in  the  main  bulkhead 
switch  panel  for  the  High  Pressure  Air  Compressor  to  the  ON 
position. 

16.  In  the  compressor  shack,  throw  the  overload  trip  arm  to  the 
ON  position  and  push  the  START  button. 

17.  After  the  compressor  has  started  runnings  place  the  unloader 
valve  in  the  RUN  position  (handle  horizontal  with  the  index 
head  pointing  towards  the  compressor.) 

18.  Adjust  the  compressor  discharge  pressure  regulator  valve  to 
give  the  desired  discharge  pressure.  (Maximum  500  psi) 

19.  Maintain  a  check  on  the  storage  tank  pressure  as  an  automatic 
unloader  for  the  compressor  is  not  provided. 

20.  Should  an  excess  amount  of  air  be  coming  out  of  the  drain  line 
from  the  trap,  it  is  probable  that  additional  water  is  required 
in  the  trap.  Water  may  be  added  when  the  compressor  is  secured 
through  the  union  connection  of  the  balance  line  until  it  is 
observed  coming  out  of  the  trap  drain. 

21.  To  secure  the  compressor  place  the  unloader  valve  in  the  START 
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position,  push  the  STOP  button  and  remove  the  oil  vent  cover* 
An  emergency  STOP  button  is  also  located  with  the  Bailey  con¬ 
troller. 

22.  Close  the  inlet  valves  to  the  storage  tanks. 

23.  Close  the  balance  line  valve. 

24.  Close  the  valve  connecting  the  separator  and  trap. 

25.  Close  the  air  supply  valve  to  the  Bailey  Controller. 

26.  Blow  down  excess  air  in  the  lines  from  the  compressor  to  the 
storage  tank  inlet  valves  by  opening  the  trap  drain  valve. 

27.  Secure  the  cooling  water. 

28.  The  drain  on  each  high  pressure  tank  and  the  inner  cooler 
should  be  blown  at  periodic  intervals  to  remove  any  accumula¬ 
tion  of  oil  or  water. 
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II.  Blow  down  procedure 

1.  Open  the  outlet  valves  on  each  high  pressure  tank. 

2.  Open  the  air  supply  valve  to  the  stagnation  pressure  Bailey 
Controller  (Supply  pressure  should  read  28  psi). 

3.  Open  the  quick  opening  valve  at  the  outlet  of  the  regulator 
va 1 ve . 

4.  Slowly  turn  the  adjusting  knob  of  the  Bailey  controller 
counter  clockwise  to  allow  a  gradual  buildup  of  pressure  in 
the  stilling  tank  to  the  desired  stagnation  pressure.  A  load¬ 
ing  pressure  of  at  least  5  psi  is  required  before  the  re¬ 
gulator  valve  will  -open. 

5.  Figure  1-A  gives  the  approximate  differential  pressure  requir¬ 
ed  to  obtain  the  maximum  Mach  number  in  the  tunnel  test  sec¬ 
tion  for  a  given  stagnation  pressure  measured  above  atmospheric 
pressure  in  inches  of  mercury. 

6.  The  time  of  each  run  may  be  prolonged  as  well  as  reducing  the 
rate  of  change  of  stagnation  temperature  by  allowing  the  com¬ 
pressor  to  continue  charging  the  storage  tanks  while  operat¬ 
ing  the  tunnel. 
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APPENDIX  B 


AIR  MOISTURE  CONTENT  CALCULATIONS. 

Assume  that  the  perfect  gas  law  relationship  applies 

PV  =  mRT 


where  P  =  gas  pressure  / bf/f 
V  =  gas  volume  f r 
m  =  lb  mass  of  gas 

R  =  universal  gas  constant  of  gas, 

T  =  gas  temperature,  °  Rankine. 

The  ratio  of  mass  water  vapor  in  the  air  to 
is  then  given  by 


f+  - ibf 
!bm  -  °R 

the  mass  of  air 


£>  if 

m  v  _  R,  T  R 

ma  ~  £  V  /?„  Pa 

Rc,  T 


where  the  subscripts  a  refer  to  air  and  v^  refers  to  the 
water  vapor.  Utilizing  the  value  of  Ra  =  53.34  and 

Rv  =  85.7 

— • '  =  0.<oZZ5  §  =  0.6225  -p' 

WQ  TOTAL  v 

Upon  exit  from  the  heat  exchanger,  the  air  is  saturated  at 
500  psia  and  60°F.  From  steam  tables  at  60  F 

Psat  r  0.2563  psia 


and 


rJL 

mQ 


t  0.6225* 


0.2563  _  -0.0003W^2?!£f?r 

500-O.25G3  “  ib.dtyo  r 
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APPENDIX  C 


FLOW  RATE  CALCULATION  PROCEDURE. 

The  mass  rate  of  flow  for  the  tunnel,  as  measured  at 
the  orifice  plate,  was  determined  by  the  procedure  out¬ 
lined  by  the  American  Society  of  Mechanical  Engineers  [5] . 
All  equations,  figures  or  table  numbers  correspond  to  num¬ 
bers  in  that  publication.  Due  to  the  removal  of  the  water 
vapor  by  the  air  drying  system,  it  is  assumed  the  air  is 
dry  for  metering  purposes.  In  addition,  due  to  the  low 
velocity  of  the  air,  the  air  is  considered  incompressible. 


For  a  representative  run: 

P  ^  =  72.08  in.  Hg. 

T  =  525°R  =  25°F 

o 

pob  =  30.84  in.  Hg. 

D  =  Pipe  diameter  =  3.120  inches, 
d  =  2.496"  =  orifice  diameter 

P  orifice  diameter  -  jJ  =  0.800 

pipe  diameter  d 

AP  =  differential  pressure  in  inches  of  H20  -  4.55 


Equation  2,  page  57,  gives  the  mass  rate  of  flow  lbm/hr. 


m  = 


359CFd2F 

a 


\!  h  r 

w 
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Where  C 


F 

F. 


h 


w 


coefficient  of  discharge 
velocity  of  approach  factor 
thermal  expansion  factor 
differential  pressure  in  inches  of  H^O 

specific  weight  of  air  at  the  inlet  side  of  the 
primary  element,  lb/ft3 


From  Figure  38,  page  67,  for  a  stainless  steel  orifice  plate 


at  65°F 

Fa  =  1.00 

From  Table  9,  page  80,  at  65°F 

?  air  =  0.075650  lb/ft3  at  29.92  in.  Hg. 
Correcting  for  increase  in  pressure 

'f  air  =  0.07565  x  =  0.0762  lb/ft3 


d2  =  (2.496  in. )2 


6.230  in. 

Assume  a  value  of  K  *  CF  =  0.80  and  calculate  an  approx¬ 

imate  m 


m  =  359  x  0.800  x  6.230  \/  4.55  x  0.0762 

■=  1053  lbm/hr. 

And  g  ^  from  Figure  35,  page  62,  for  65°F  air  =  .0000122 
Calculate  Reynolds  number  from  Equation  8. 


«D 


0.004244  m  _  0. 004244  x  1053  =.  117,600 


Dg 


A 


3.120  x  0.0000122 


From  Table  5,  page  26,  or  from  Figure  1-C  of  this  Appendix 
for  =  0.800,  and  R  =  117,600  the  flow  coefficient 

is  K  =  0.8062 

Corrected  calculation  of  flow  rate,  therefore  is 


m  =  359~x--0. 8062  x  6.230  *V  4.55  x  0.0762 

_  1060  lbm/hr.  =  0.2945  lbm/sec. 
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The  ideal  mass  flow  rate  at  the  nozzle  throat  is  determined 
by  Fliegner's  formula  [2]. 

A*P 


m; 


Where 


M  = 


1.00 


0.532 


^  T 1 

r  o 

=  ideal  flow  rate  Ibm/sec. 

=  tunnel  throat  area  in  square  inches  with 


=  stagnation  pressure  psia 
=  stagnation  temperature  in  °Rankine 


0.532  x  0.367  x  72.08  x  0.491 

-\f~525~ 


=  0.301  lbm/sec. 


The  nozzle  discharge  coefficient  CD  is  defined  as  the 
ratio  of  the  actual  nozzle  flow  to  the  flow  calculated 
from  the  isentropic  laws  for  the  initial  and  final  pres¬ 
sures  of  the  actual  nozzle.  Based  on  this  definition  the 
coefficient  of  discharge  is 


m 


D 


actual 


m 


ideal 


0. 2945 

0.  301 


0.969 


Figure  1-C  presents  the  Flow  Coefficient  K  plotted 
vs.  the  Reynolds  Number  for  a  flat  plate  orifice  with 
pressure  taps  at  D  and  1/2  D  for  a  three  inch  pipe.  This 
figure  presents  the  data  of  Table  5  (5)  graphically  for 
values  of  d/D  equal  to  0.580  and  0.800,  Figure  25  presents 
graphically  the  nozzle  discharge  coefficient  plotted  vs, 
the  stilling  tank  stagnation  pressure.  Values  of  were 
determined  by  the  use  of  Fliegner's  formula  and  by  the 
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procedure  outlined  by  the  American  Society  of  Mechanical 
Engineers  (5)  as  used  in  this  Appendix. 
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APPENDIX  D 


DESIGN  OF  DIFFERENTIAL  PRESSURE  TRANSDUCER. 

Differential  pressure  transducers  were  originally  to  have  measured 
the  differential  pressure  between  the  tunnel  pressure  and  the  stagna¬ 
tion  pressure.  This  method  of  measurement  was  discarded  due  to  the 
difficulties  of  fabricating  similar  output  transducers. 

A  successful  transducer,  Figures  2-D  and  3-D,  was  made  based  on  the 
principle  of  the  deflection  of  a  circular  plate  held  firmly  at  the  edge 
>]  • «  was  desired  that  each  psi  of  differential  pressure  produce  at 
least  10  microinches  per  inch  of  strain.  Transforming  the  equation  of 
stress  at  the  center  of  the  clamped  plate  to  the  equation  of  strain  (end 
of  this  appendix)  and  substituting  the  appropriate  values  the  plate  thick- 
ness  was  calculated  as  being  0.0336  inches  times  the  radius  of  the  plate. 
One  A-7  strain  gage  was  located  at  the  center  of  each  side  of  the  plate 
and  the  cumulative  output  of  the  two  gages  was  read  on  a  Baldwin  strain 
indicator.  This  cumulative  output  was  thus  expected  to  be  approximately 
20  microinches  per  inch  of  strain  per  psi  of  differential  pressure. 

The  plate  thickness  actually  used  was  0.032  inches  due  to  the 
availability  of  this  uniform  thickness  sheet  stock.  This  plate  was 
welded  between  two  hollow  cylinderical  shells  of  1.5  inches  inside 
diameter  and  2.25  inches  outside  diameter.  While  this  combination  of 
plate  thickness  and  effective  plate  diameter  exceeded  the  calculated 
plate  thickness  to  plate  diameter  ratio,  the  combined  output  still  ex¬ 
ceeded  the  desired  minimum  of  10  microinches.  The  open  ends  of  the 
cylinderical  shells  were  threaded  to  receive  a  plug  which  contained  a 
pressure  tap.  An  0-ring  seal  was  used  between  the  plug  and  the  shell 
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of  the  transducer. 


One  leg  of  each  strain  gage  was  led  through  the  shell  wall  and 
the  drilled  hole  secured  with  epoxy  cement .  The  other  leg  of  the  gage 
was  soldered  to  the  inside  of  the  shell  which  served  as  a  common  ground 
between  the  two  gages.  The  two  protruding  wires  from  the  gages  and  a 
common  ground  wire,  soldered  to  the  outside  of  the  shell,  were  then 
connected  to  the  Baldwin  strain  gage  indicator. 

The  completed  transducer  was  calibrated  on  a  hydraulic  dead  weight 
tester.  One  end  of  the  transducer  was  vented  to  the  atmosphere  and  the 
other  cavity  filled  with  oil  and  connected  to  the  dead  weight  tester. 

A  pressure  differential  ranging  from  0  to  100  psi,  in  increments  of 
10  pounds^  was  then  imposed  with  strains  recorded  while  loading  and  un¬ 
loading  the  tester. 

A  linear  plot  of  pressure  differential  vs.  strain  was  obtained  for 
several  transducers  but  on  others,  the  welding  caused  a  slight  permanent 
buckle  in  the  plate  and  the  strain  indications  were  not  reproduceable 
for  a  given  pressure  differential.  Due  to  the  time  involved  in  making 
and  testing  each  transducer,  it  was  felt  that  too  much  unproductive  time 
would  be  spent  on  making  transducers  that  could  not  be  used;  so  this 

■  t  • 

means  of  pressure  measurement  was  abandoned  in  favor  of  manometer 
pressure  measurements. 

One  test  transducer  w^s  gltled  together  with  Eastman  910  glue  to 
eliminate  the  problems  caused  by  the  buckling  of  the  diaphragm  due  to 
welding.  A  centering  hole  and  shoulder  was  made  in  the  cylinderical 
shells,  the  diaphragm  was  then  glued  to  both  shells  under  light  pressure. 
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The  strain  gages  were  attached  to  each  side  of  the  diaphragm  prior  to 
joining  the  two  end  shells.  This  transducer  gave  a  linear  plot  of  18,52 
microinches  per  inch  of  strain/psi  of  pressure  differential  when  tested 
on  the  hydraulic  dead  weight  tester  but  it  did  not  respond  to  less  than 
3/4  of  an  inch  of  mercury  pressure  differential  when  checked  against  a 
manometer  reading  of  the  same  pressure  differential  with  the  tunnel  in 
operation,. 

The  stress  distribution  in  a  uniformly  loaded  diaphragm  with  clamp¬ 
ed  edges  is  given  in  Figure  1-D  below  M  . 


Figure  1-D  Stress  Distributions  of  Diaphragm 


The  tangential  stress  at  the  center  of  the  plate  is  given  by  the 


equation 


S  = 


Where 


3U±£]_P£ 

8  i. 2 


S  =  stress  in  psi. 

=  Poisson  ratio. 

jO  =  Pressure  differential  in  psi. 

P  =  Radius  of  diaphragm  in  inches. 

■£  =  Diaphragm  thickness  in  inches. 

The  corresponding  value  of  strain  is  obtained  from 


(1) 


6  -  -  -  f  4 

^  a  *  e 


(2) 


where  £  is  the  strain  and  £  =  modulus  of  rigidity  of  diaphragm 
material.  Substituting  equation  (2)  into  equation  (1) 


€  - 


3P  r*  (I —  H 

BE  tz 


(3) 


Utilizing  the  required  minimum  value  of  £  =  10  7^“  /Ps^ 

differential,  =  0.3  and  E  =  30  x  10^  psi  for  the  values  of 

steel  the  ratio  of  the  diaphragm  thickness  to  its  radius  may  be  obtained. 


to 


-5 


3'  h  tjU  /-a  3l  ) 

8  *  3*  ‘o'1  *  lz 


i  -  O.o  33(o 


w 


For  a  diaphragm  radius  of  0.75  inches  a  plate  thickness  of  0.0252 


inches  is  required. 
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Figure  2-D.  Exterior  View  of  Welded  Pressure  Transducer. 


Figure  3-D.  Interior  View  of  Welded  Pressure  Transducer. 
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APPENDIX  E 


GENERAL  THEORY  OF  THE  SCHLIEREN  METHOD  OF  OPTICAL  INVESTIGATION  OF 
COMPRESSIBLE  FLUID  FLOW  [2]  . 

Apart  from  the  conventional  methods  of  experimentally  investigating 
flow  patterns  oy  means  of  pressure  and  velocity  surveys,  compressible 
gas  flows  lend  themselves  particularly  well  to  optical  methods  of  in¬ 
vestigation. 

Fundamentally,  the  optical  methods  in  common  use  (the  Interfero¬ 
meter,  the  Schlieren,  and  the  shadowgraph)  depend  on  one  of  the  two 
physical  phenomena:  (1)  the  speed  of  light  depends  on  the  index  of 
refraction  of  the  medium  through  which  it  passes,  and  the  index  of  re¬ 
fraction  of  a  gas  in  turn  depends  upon  its  density;  and,  as  a  conse¬ 
quence  of  this  first  phenomenon,  (2)  light  passing  through  a  density 
gradient  in  a  gas  (and  therefore  through  a  gradient  in  index  of  re¬ 
fraction)  is  deflected  in  the  same  manner  as  though  it  were  passing 
through  a  prism.  In  a  high-speed  gas  flow  the  density  changes  are 
sufficiently  large  to  make  these  phenomena  sizeable  enough  for  optical 
observation. 

The  interferometer,  based  on  phenomenon  (1),  measures  directly 
changes  of  density,  and  is  primarily  suited  for  quantitative  determina¬ 
tion  of  the  density  field. 

The  Schlieren  method,  based  on  phenomenon  (2),  measures  density 
gradients.  Although  it  is  theoretically  adaptable  to  quantitative  use, 
it  is  inferior  to  the  interferometer  in  this  respect,  and  its  greatest 
utility  is  in  giving  an  easily  interpretable  picture  of  the  flow  field 
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together  with  a  rough  picture  of  the  density  variations  in  the  flow. 

The  shadowgraph  method,  also  based  on  phenomenon  (2),  measures 

the  second  derivative  of  the  density  (i.e.,  the  first  derivative  of 

the  density  gradient).  Therefore  it  makes  visible  only  those  parts 

of  the  flow  where  the  density,  gradients  change  very  rapidly,  and  it 

has  found  its  greatest  utility  in  the  study  of  shock  waves. 

Of  the  three  methods  mentioned,  the  interferometer  yields  the 

most  information,  and  the  shadowgraph  the  least.  On  the  other  hand, 

% 

the  interferometer  is  the  most  costly  and  the  most  difficult  to  operate^ 
whereas  the  shadowgraph  is  the  least  costly  and  the  easiest  to  operate. 
Each  method  therefore  has  its  own  useful  niche  in  experimental  workj, 
and  the  choice  of  method  depends  on  the  nature  of  the  investigation. 
Figure  1-E  below  shows  one  of  the  numerous  type  of  Schlieren 
(striae)  arrangements.  Light  from  a. uniformly  illuminated  line  source 


rH 


Figure  1-E  Schematic  Schlieren  System 
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of  small  but  finite  width  is  collimated  by  the  first  len9  and  then 
passes  through  the  test,  section.  It  is  then  brought  to  a  focus  by 
the  second  lens  and  projected  on  the  screen.  At  the  focal  point, 
where  there  exists  an  image  of  the  source,  there  is  introduced  a  knife- 
edge  which  cuts  off  part  of  the  light.  With  no  flow  in  the  test  sec¬ 
tion  the  knife-edge  is  usually  adjusted  so  as  to  intercept  about  half 
the  light,  and  the  screen  is  uniformly  illuminated  by  the  portion  of 
the  light  escaping  the  knife-edge.  When  the  flow  is  established  in 
the  test  section  assumed ,  for  simplicity,  to  be  two  imensional,  with 
oach  light  ray  passing  through  a  path  of  constant  air  density  ar.y  light 
i ay  passing  through  a  region  in  which  there  is  a  density  gradient  nor¬ 
mal  to  the  light  direction  will  be  deflected  as  though  it  has  passed 
through  a  prism.  Therefore,  depending  on  the  orientation  of  the  knife- 
edge  with  respect  to  the  density  gradient,  and  on  the  sign  of  the  den¬ 
sity  gradient,  more  or  less  of  the  light  passing  through  each  part  of 
the  test  section  will  escape  the  knife-edge  and  illuminate  the  screen. 

Thus  the  S:hlieren  system  makes  density  gradients  visible  in  terms 
of  intensity  of  illumination.  A  photographic  plate  at  the  viewing 
screen  records  density  gradients  in  the  test  section  as  different  shades 
of  gray. 

For  the  Schlieren  figure  drawn  it  is  assumed  that  the  flow  in  the 
test  section  is  parallel  and  in  the  x,  y-plane,  and  that  the  light 
passes  through  the  test  section  in  the  z-direction. 

Since  the  speed  of  a  wave  front  of  light  varies  inversely  with 
the  index  of  refraction  n  of  the  medium  through  which  the  light  travels, 
it  follows  that  a  given  wave  front  will  rotate  as  it  passes  through  a 
gradient  in  n,  Accordingly,  the  normal  to  the  wave  front  will  follow 
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a  curved  path;  since  this  normal  is  what  we  mean  by  the 
light  ray,  the  latter  is  refracted  as  it  passes  through  the 
density  gradient.  Noting  that  n  is  nearly  unity  for  gases, 
it  may  be  shown  that 
I 

—  r  Gradient  n 
R  * 

where  R  is  the  radius  of  curvature  of  the  light  ray.  The 
total  angular  deflection  of  the  ray  in  passing  through 

the  test  section  of  width  L  is  therefore  given  by 

*7“--  “  L  g  ra  d  n 

Resolving  this  into  Cartesian  components,  and  connecting 
the  index  of  refraction  with  the  gas  density  through  the 
empirical  Gladstone-Dale  equation  which  states  that 


M=!  -  /< 

9  g-d 

where  K_  _,  the  Gladstone-Dale  constant,  is  constant  for  a 
G 

given  gas,  we  get  the  angulax  deflections  of  the  light  in 
the  x-  and  y-directions 


L  K 


jj? 


If  the  knife-edge  is  aligned  in  the  y-direction  (i.e,, 
normal  to  the  flow)  as  in  Figure  1-E,  only  deflection  in 
the  x-direction  will  influence  the  light  passing  the  knife- 
edge.  Hence  density  gradients  in  the  x-direction  will  be 
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made  visible,  but  gradients  in  the  y-direction  will  not  be 
visible.  Similarly,  if  the  knife-edge  is  parallel  to  x, 
only  gradients  in  the  y-direction  will  become  visible. 

In  interpreting  a  Schlieren  photograph  of  a  two  dimen¬ 
sional  flow,  it  is  convenient  to  imagine  that  the  photograph 
represents  a  view  of  a  relief  map  of  the  flow  field  in  which 
the  vertical  elevation  is  the  gas  density.  If  the  map  is 
illuminated  by  glancing  side  light  nearly  parallel  to  x,  the 
pattern  of  shadows  and  brightness  will  correspond  to  a 
Schlieren  photograph  with  the  knife-edge  in  the  y-direction. 
Similarly,  illumination  in  the  y-direction  corresponds  to 
placing  the  Schlieren  knife-edge  parallel  to  the  x-direction. 

Figures  14b  and  14c  are  Schlieren  photographs  of  the 
installed  test  section  and  show  the  two  different  effects 
obtainable  with  orientations  of  the  knife-edge.  In  Figure 
14b  (knife-edge  horizontal)  light  areas  on  the  upper  half 
correspond  to  dark  areas  on  the  lower  half,  and  vice  versa; 
the  boundary  layers  on  the  walls  are  clearly  visible,  since 
with  this  knife-edge  orientation  the  density  gradients  nor— 
mal  to  the  flow  are  made  visible.  In  Figure  14c  (knife  edge 
vertical)  the  upper  and  lower  halves  look  alike;  the  bound¬ 
ary  layers  are  not  visible,  except  aft  of  the  shock  wave 
which  is  inneracting  with  the  boundary  layer. 

The  quality  of  focus  of  the  Schlieren  pictures  is  be¬ 
lieved  to  be  poor  due  in  part  to  the  present  arrangement  of 
the  Schlieren  system  installed  for  the  reasons  as  outlined 
below  by  Holder  and  North  [7] : 

Although  it  is  usually  desirable  to  arrange  for  the 
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beam  of  light  passing  through  the  working  section  to  be 
parallel,  nonparallel  beams  are  occasionally  used  and  an 
arrangement  for  producing  this  type  is  sketched  in  Figure 
2-Eo  Only  one  mirror  or  lens  is  required  (apart  from  the 
focussing  lens) ,  but  this  advantage  is  to  some  extent  off¬ 
set  because  its  diameter  must  be  larger  in  relation  to  the 
field  to  be  examined  than  if  a  parallel  beam  is  used. 


Figure  2-E.  Coincidence  Schlieren  System  using 
a  Plane  Mirror  or  knife-edge 


In  addition  to  employing  a  nonparallel  light  beam, 
the  mirror  system  involves  a  double  passage  of  the  beam 
across  the  working  section.  The  source  is  placed  at  the 
center  of  curvature  of  the  mirror ,  and  part  of  the  image 
of  the  source  is  reflected  on  to  the  viewing  screen  by  a 
small  plane  mirror  held  near  the  source.  Because  the 


64 


source  and  its  image  are  very  close  to  each  other,  the 
arrangement  is  often  referred  to  as  the  "coincidence"  system. 
If  the  image  of  the  source  is  deflected  in  a  direction  per¬ 
pendicular  to  the  edge  of  the  plane  mirror,  the  illumination 
on  the  viewing  screen  changes;  the  edge  of  the  mirror  thus 
plays  the  same  part  as  the  knife  edge  in  a  Toepler  appa¬ 
ratus. 

Since  the  source  and  its  image  are  not  truly  coincident, 
the  paths  of  the  incident  and  reflected  rays  across  the 
working  section  differ  even  in  the  absence  of  deflections 
produced  by  density  gradients.  This  has  an  adverse  effect 
on  the  image  quality  (in  addition  to  any  effect  associated 
with  the  use  of  non-parallel  light)  since  two  slight  offset 
images  are  produced.  (See  Figure  16c) 

These  images  may  be  made  to  coincide  at  the  expense  of 
considerably  reduced  illumination  by  using  the  arrangement 
sketched  in  Figure  3-E. 


Figure  3-E.  Coincidence  Schlieren  System  using 
Graded  Filter 
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Here  a  graded  filter  is  placed  at  45  degrees  to  the  axis  at 
the  center  of  curvature  of  the  mirror,,  A  condenser  lens  is 
used  to  produce  an  image  of  the  source  on  the  surface  of 
the  filter ,  and  part  of  the  light  from  this  image  is  re¬ 
flected  into  the  mirror.  The  image  of  the  source  formed  by 
reflection  from  the  mirror  is  made  to  coincide  with  the 
original  image  on  the  surface  of  the  filter,  and  the  light 
transmitted  by  the  filter  is  arranged  to  fall  on  the  view- 
ing  screen  in  the  usual  way.  The  source  and  its  image  may 
thus  be  made  exactly  coincident  so  that,  in  the  absence  of 
density  gradients,  the  incident  and  reflected  beams  coin¬ 
cide.  This  will  no  longer  be  the  case  when  the  rays  are 
deflected  by  density  gradients,  and  some  loss  of  image 
sharpness,  or  difficulty  in  interpreting  the  observations 
may  then  result. 

A  further  difficulty  arises  because  only  one  of  the 
two  images  formed  by  the  incident  and  the  reflected  light 
can  be  in  focus.  It  is  usually  best  to  focus  the  image 
associated  with  the  reflected  light  since  the  incident  rays 
are  subject  to  further  refraction  during  their  second  pass¬ 
age  across  the  wind  tunnel. 

The  principal  advantage  of  the  system  is  that  the 
double  passage  of  the  light  beam,  and  the  fact  that  the 
source  image  is  formed  at  the  center  of  curvature  of  the 
mirror  instead  of  its  focus,  give  increased  sensitivity. 

For  example,  in  a  two-dimensional  flow  containing  a 
region  of  density  gradient  of  sufficient  extent  for  the 
incident  and  reflected  ray  to  remain  within  it,  the. 
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angular  deflection  of  the  ray  will  be  doubled  by  the  use  of 
such  a  system.  Moreover,  the  displacement  of  the  image  of 
the  source  for  a  given  angular  deflection  will  be  twice  as 
large  when  the  image  is  at  the  center  of  curvature  as  when 
it  is  at  the  focus.  Coincidence  systems  may  thus  be  val¬ 
uable  where  extreme  sensitivity  is  required  (for  example, 
tunnels  working  with  very  low  density) ,  although  for  or¬ 
dinary  use  they  are  undesirable  for  the  reasons  mentioned 
above. 

Part  of  the  increased  sensitivity  associated  with  a 
double-passage  system  may  be  retained,  without  the  diffi¬ 
culties  associated  with  the  use  of  non  parallel  light  by 
means  of  the  arrangement  sketched  in  Figure  4-E.  Here  a 
parallel  beam  of  light  is  produced  by  a  lens  or  mirror  and 
after  crossing  the  working  section  it  is  reflected  back  by 
means  of  a  plane  mirror. 


Figure  4-E.  Schlieren  System  with  Parallel  Light  Beams 
Crossing  Tunnel  Twice  with  Increased  Sensi¬ 
tivity. 
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Tinsley  Laboratories  Inc.  furnished  the  basic  components  for 
the  coincidence  Schlieren  system  as  well  as  the  design  as 
shown  in  Figure  11.  One  unit  was  the  six  inch  spherical 
mirror  with  a  focal  length  of  80  inches  mounted  on  a  small 
pedestal.  The  other  unit,  mounted  on  a  wooden  frame  with 
leveling  jacks,  consisted  of  the  light  source,  knife  edge, 
surface  mirror  and  encased  in  a  light  tight  wooden  box  the 
camera  lens,  folding  surface  mirrors  and  a  six  inch  ground 
glass  view  plate.  These  were  mounted  on  a  portable  stand 
as  shown  in  Figure  9. 

An  Ilex  Optical  Company  Number  3  Acme  Synchro  Shutter 
was  installed  between  the  knife  edge  and  the  camera  lens  by 
means  of  a  machined  adapter  fitted  over  the  lens  mounting 
with  the  shutter  bolted  to  it.  This  arrangement  allowed 
all  of  the  light  passing  over  the  knife  edge  to  enter  the 
camera  lens  system  but  otherwise  maintained  the  light  tight 
integrity  of  the  system. 

The  six  inch  ground  glass  view  plate  was  removed  and 
the  back  of  an  Ilex  Optical  Company  Korona  camera  consist¬ 
ing  of  a  four  inch  by  five  inch  ground  glass  view  plate 
and  film  pack  holder  fitted  in  its  place. 

The  procedure  to  align  the  Schlieren  system  is  as 
follows  referring  to  Figure  11: 

1.  Locate  the  spherical  mirror  as  close  to  the  test 
section  window  as  practical. 

2.  Locate  camera,  light  source,  and  knife-edge 
assembly  such  that  the  light  source  lies  on  an  axis 
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perpendicular  to  the  window.  The  light  source  should  be 
approximately  80  inches  from  the  spherical  mirror  surface . 
This  may  be  done  by  moving  the  spherical  mirror  in  the  slots 
provided  on  the  stand. 

3.  With  the  light  source  "turned-On"  adjust  the 
spherical  mirror  to  image  the  light  source  at  the  knife- 
edge.  If  the  image  is  not  in  focus  with  the  knife-edge 
assembly  at  approximately  the  middle  of  its  slide  travel 
move  the  spherical  mirror  forward  or  back  to  correct  the 
image  focus.  (A  small  black  card  is  helpful  in  locating 
images  and  the  light  beam  directions. ) 

4.  Adjust  the  knife-edge  screw  such  that  the  "Knife" 
does  not  intercept  the  light  source  image.  The  light  should 
then  pass  through  the  camera  lens,  and  an  image  of  the  test 
region  should  appear  on  the  view  screen.  If  the  view  screen 
image  is  not  centered,  adjust  the  small  1st  surface  mirror 
and  the  spherical  mirror  such  that  the  light  passes  through 
the  center  of  the  camera  lens  and  aperture  of  the  knife- 
edge  assembly.  If  the  view  screen  image  is  yet  not  centered, 
adjustment  of  the  three  small  screw  heads  at  the  back  of  the 
camera  will  provide  the  final  centering  adjustment. 

5.  The  knife-edge  can  now  be  freely  adjusted  to  provide 
the  required  Schlieren  sensitivity. 

NOTE:  All  mirrors  are  front  surface  aluminized  and 

should  not  be  touched  on  the  reflecting  surface. 
They  should  be  cleaned  by  dusting  with  a  qual¬ 
ity  camels  hair  brush. 
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The  best  photographs  of  the  flow  field  were  obtained 
using  the  following  camera  settings: 

1.  Shutter  speed  set  at  1/200  of  a  second. 

2.  Lens  opening  at  f-11.  This  value  has  nothing  to 
do  with  the  actual  exposure  time  of  the  film. 

This  value  of  the  lens  opening  allows  all  of  the 
light  beam  to  enter  the  camera  and  helps  reduce 
stray  light  entering  the  lens  system  and  fogging 
the  film. 

3.  Royal  Pan  film  with  a  tungsten  film  speed  rating 
of  160  was  used. 

4.  Developing  time  of  eight  minutes  in  total  darkness 

in  a  1  :  1. dilute  mixture  of  Eastman  DK  50  devel- 

o 

oper  and  water  at  68  F. 

Colored  Schlieren  views  of  the  flow  field  may  be  ob¬ 
tained  by  replacing  the  knife-edge  assembly  with  the  color 
grid  purchased  from  Tinsley  Laboratories.  Colored  pictures 
were  not  attempted  as  a  means  of  printing  and  developing 
the  color  film  was  not  available. 

The  laght  source  may  be  either  a  General  Electric 
Company  bulb  number  1968  or  number  1924,  or  a  Westinghouse 
bulb  number  1633.  All  three  bulbs  produce  a  line  source 
of  light.  If  the  General  Electric  bulbs  are  used,  a  dif¬ 
ferent  mount  must  be  made  for  the  bulb. 

Both  the  light  source  and  the  knife  edge  may  be  rotated 
from  the  horizontal  position  to  the  vertical  position  to 
obtain  the  desired  direction  of  the  density  gradient  of  the 
flow  field. 
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For  the  reasons  already  outlined,  the  focus  of  this 
Schlieren  system  is  not  as  sharp  as  could  be  obtained  with 
a  two  mirror  one  pass  light  beam  system, 

A  metal  frame  was  fitted  over  one  of  the  Schlieren 
windows  and  small  copper  wires  were  soldered  into  position 
at  one  inch  intervals,  starting  at  the  throat  position,  to 
give  horizontal  dimensions. 
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APPENDIX  F 


SAMPLE  CALCULATIONS  FOR  EXPERIMENTAL  DATA  AND  THEORETICAL 
PREDICTIONS. 

PREDICTION  OF  ADIABATIC  MACH  NUMBERS  FOR  NOZZLE  AND  CON¬ 
STANT  AREA  TEST  SECTION  OPERATING  WITHOUT  SHOCK. 

Flow  is  assumed  to  be  isentropic  to  the  nozzle  throat 
and  adiabatic  thereafter.  For  ease  of  calculations,  the 
nozzle  contour  as  given  in  Figure  19a  was  assumed  to  con¬ 
sist  of  five  straight  line  segments  as  listed  in  Table  X. 

Utilizing  the  influence  coefficients  of  table  8.2  (2) 
multiplied  by  M  and  designating  the  product  F  as  outlined 
by  Shapiro  (2)  for  the  combined  effects  of  area  change  and 
friction  results  in  the  following  equation 

dM2  =  F  _|_  F^4f  dx  . 

A  A  *  f  D 

The  technique  of  isoclines  as  given  in  Chapter  8  of  Shapiro 
(2)  was  the  means  of  arriving  at  a  solution  to  this  differ¬ 
ential  equation. 

The  parameters  of  the  area  A  and  the  hydraulic  dia¬ 
meter  D  may  be  expressed  in  terms  of  the  tunnel  dimensions. 
Area  of  the  tunnel  at  any  point  is  the  product  of  the  tunnel 
width  w  and  the  tunnel  height  h  at  the  point  in  question. 

As  the  tunnel  is  of  constant  width,  the  area  at  any  point  is 
proportional  to  the  height  at  that  point  and  the  derivative 
of  the  area  proportional  to  the  derivative  of  the  tunnel 
height.  The  hydraulic  diameter  is  defined  as  4  times  the 
area  divided  by  the  wetted  perimeter  and  expressed  in 
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tunnel  dimensions  is 


4  w  h 
2  (w+  h) 


D  =  4  x  area 

wetted  perimeter 


2  w  h 

w+  IT 


As  the  tunnel  width  w  is  constant  throughout  the  length  of 
the  tunnel  at  0.500  inches  the  hydraulic  diameter  may  be 
expressed 


D 


o  (0. 500)  h 

1  o:for+h 


2  h 
lt2h 


As  assumed,  the  nozzle  segment  heights  will  vary  lin¬ 
early  in  accordance  with  the  equation  for  a  straight  line 
h  =  mx  +  hQ  where  m  is  the  slope  of  the  nozzle  segment 
at  some  distance  x  downstream  of  the  nozzle  throat  and  the 
nozzle  throat  having  an  initial  height  of  hQ. 


at  x 


h  =  mx+  h 

•  o 

dh  ~  mdx  ;  dx  =  — 

m 

L,  h  =  h^  so  the  slope  m 


*>  L 


Combining  these  relationships  and  substituting  them  into 
the  original  differential  equation 


dM2  - 


A  TT  +  Ff(4f)  2 


a  -f  \  L  (  w  •+  h , 


w  h  (h^  -  h  ) 


dh 


dM2  = 


dM 


F  F  I4f)  L 

A  ,  f 

YT  + 


2  w  (hx  -  h  ) 


dh 


F  F, (4f ) L 

-A  +  _ _ 

h  2  w  (h,  -  h 


1  ^  )  dh 

h 


(i+“> 


The  equation  is  now  in  the  form  so  that  by  suitable  numer¬ 


ical  values  of  h  and  M  being  substituted  \n:he  slope  of  the 


,Q;V 
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Mach  number  squared  is  obtained.  Values  of  are  assumed 

in  addition  to  values  of  M  that  are  thought  to  bracket  the 

2 

actual  Mach  number.  On  a  plot  of  M  vs,  hi  the  calculated 

2 

values  of  the  slope  dM  /dh  are  plotted  at  the  coordinate 

points  of  each  assumed  value  of  M  and  h.  Sufficiently  small 

increments  of  h  must  be  utilized  to  provide  the  necessary 

isoclinic  field,  A  smooth  curve,  starting  tangent  to  be 

2 

the  value  of  dM  /dh  at  the  known  value  of  h^  and  its  corres- 

o 

ponding  value  of  M  is  then  drawn  through  the  field  of  iso¬ 
clines  tangent  to  a  reference  a socline  at  each  value  of  h 
the  curve  passes  through.  In  addition  to  the  physical  data 
listed  in  the  previous  table  a  value  of  the  average  friction 
factor  (4f)  equaling  0,01  was  assumed  with  an  initial  Mach 

number  of  1,00  existing  in  the  nozzle  throat. 

2 

A  value  of  M  is  determined  from  the  junction  of  the 
smooth  curve  and  the  height  at  the  segment  end  point.  This 
calculated  Mach  number  is  then  used  as  the  initial  Mach 
number  for  the  second  nozzle  segment.  Continuation  of  this 
process  resulted  in  the  adiabatic  Mach  number  at  the  noz¬ 
zle  exit.  Sample  calculations  for  segment  2  are  given  in 
tabular  form  in  Table  XII  and  are  plotted  on  figure  20, 
Straightforward  use  of  Table  42  ( 3)  for  Fanno  line 
procedure  was  used  to  predict  the  resultant  Mach  number  at 
the  end  of  the  constant  area  section  of  the  tunnel.  Iso¬ 
clinic  procedures  could  have  been  used  but  this  method 
was  discarded  as  being  time  consuming  compared  to  the  use 
of  the  tables.  The  following  information  was  assumed  or 
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utilized  for  determining  the  exit  Mach  number:  of  f 


T 


4f 


L 


D 

Mi 


0,01  and  constant  in  value  throughout  the 
section 

5.75  inches,  the  tunnel  section  length 
0.775  inches,  the  hydraulic  diameter 
Mach  number  at  the  nozzle  exit  equal  2.17 
Mach  number  at  the  section  exit 


Entering  Table  42(3)  with  Mi  of  2.17,  a  value  of  4fL  _  /D 

-L  IHaX 

was  read  as  0.35298  for  use  in  the  following  equation: 


0.35298  *  (0.01)  (5.75)  _  0.27878 

0.775 


Table  42(3)  was  then  entered  using  this  value  of  4fL/D  and 
the  Mach  number  at  the  end  of  the  test  section  determined 
as  1.914. 


PREDICTION  OF  ADIABATIC  MACH  NUMBER  WITH  A  NORMAL  SHOCK 
STANDING  IN  THE  TUNNEL  NOZZLE. 

Flow  was  assumed  to  be  adiabatic  in  the  nozzle  until 
the  normal  shock  occurred  at  a  point  midway  between  pres¬ 
sure  taps  4  and  5.  M^  the  Mach  number  upstream  of  the 
normal  shock  was  determined  as  2.0  as  read  from  figure  20, 
which  is  a  plot  of  the  adiabatic  Mach  numbers  as  calculated 
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vs.  tunnel  length  based  on  previous  calculations  in  this 
appendix.  Entering  Table  48(3)  with  this  value  of  M  ,  the 

X 

Mach  number  M^.  downstream  of  the  normal  shock  was  obtained 
as  0.57735.  M^.  thus  obtained  was  used  as  the  initial  value 

of  M  for  the  isocline  technique  for  determining  Mach  numbers 
for  points  downstream  of  the  normal  shock  in  the  nozzle  and 
in  the  constant  area  section  of  the  tunnel.  The  nozzle  area 
and  the  corresponding  nozzle  height  were  read  from  figure 
19a  at  the  shock  location  assumed  to  be  between  pressure 
taps  4  and  5.  Nozzle  height  at  the  assumed  shock  location 
is  1.30  inches. 

As  before,  the  remaining  portion  downstream  of  the  shock 
in  the  diverging  nozzle  was  assumed  to  consist  of  the  same 
straight  line  segments.  The  initial  section  was  taken  as 
the  remaining  portion  of  the  original  segment  2  and  was 
designated  as  segment  2a.  All  previous  assumptions  con¬ 
cerning  the  average  friction  factor  and  the  physical  data 
remained  constant  as  outlined  in  Table  XI. 

A  nozzle  exit  Mach  number  was  found  to  be  0.429.  The 
relationship  of  4fL/D  at  the  inlet  and  exit  of  the  constant 
area  section  was  utilized  to  determine  the  Mach  number  at 
the  end  of  the  constant  area  section  as  was  done  in  pre¬ 
vious  calculations  in  this  appendix.  This  exit  Mach  number 
was  found  to  be  M  =  0.434. 

Theoretical  values  of  P/P  .  at  any  point  downstream  of 

ol 

the  shock  are  obtained  by  multiplying  the  value  of  P/P0 
corresponding  to  the  Mach  number  at  that  point  by  the  value 
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of  P  /P  obtained  from  Table  48  (3)  for  the  Mach  number 

at  which  the  shock  occurred.  Based  on  the  data  in  Table 

IX,  it  was  assumed  the  shock  occurred  at  a  Mach  number  of 

2.000  with  P  /P  having  a  value  of  6.7208.  Adiabatic 
oy  ox 

Mach  numbers  and  pressure  ratios  obtained  by  this  procedure 
are  presented  graphically  in  Figures  22  and  23. 

Table  XVI  is  a  summary  of  the  adiabatic  Mach  numbers 
predicted  for  the  tunnel  operating  with  and  without  a  nor¬ 
mal  shock  standing  in  the  tunnel  nozzle  with  a  constant 
area  section. 

TABLE  XVI  -  Summary  of  Adiabatic  Mach  Numbers 
in  Nozzle 

Pressure  tap  Adiabatic  Mach  Adiabatic  Mach 

Number  Number 

No  Shock  Shock  in  Nozzle 


2 

1.355 

1.355 

5 

2.070 

0.529 

7 

2.180 

0.432 

9 

2.195 

0.429 

12 

2.170 

0.429 

32 

1.914 

0.434 

(end  of  constant 
area  section) 
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DETERMINATION  OF  DIFFUSER  HEIGHT  TO  PRODUCE  CONSTANT  MACH 
NUMBER  IN  TEST  SECTION  INCLUDING  EFFECT  OF  FRICTION* 

Assume  that  the  diffuser  height  will  increase  linearly* 
Point  t  is  taken  as  the  nozzle  throat;  point  2,  the  end  of 
the  nozzle;  and  point  3  the  maximum  diffuser  height  loca¬ 
tion. 

% 

Utilizing  the  influence  coefficients  for  constant 
specific  heat  and  molecular  weight  from  Shapiro's  Table 
8.2  [2] 

^ w*) Sd*  -  q 

~p]z  "  j-M8-  '  A  /-A72  D 


Area  of  tunnel  -  width  x  height  =  w  h  and  as  tunnel  is  of 
constant  width 


dA  Z  dh  qnd  ^ 


h 


The  hydraulic  diameter  D  = 


4  (area) 

Wetted  perimeter 


4  w  h 

2  (w  +  h) 


2  w  h.  Substituting  the  values 
w  -f  h 


pjJA  -  p  dA  ~  4££j±d*,(0J+h ) 
c  /\  c  h  2o)h  x 
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Combining  terms 


4  wdh_ 

oj  +  h 


=  4 -fjk  Msd'A 


Integrating  between  designated  limits 


h  — 

cu/n  (u)+h)\ ’  =  Jk  M*-f  (X3~ X*) 

hZ 

(n  /  OJ-b  ha  )  M 

11  1  LU  +  h  %  /  "  ^ 


clearing  and  solving  for 

■kMaf(X*,~Xz) 

/L,  =  —  uj  -j-  ( uJ+  h )  & 


Substituting  values  of  w  =  0.5  inches;  -  1.718  inch; 

M2  =  2.17;  X3  -  X2  =  5.33  inches,  Figure  20. 

Xf(a  /7)2(a^g5X433) 

h3  =  0.5  +  (2.218) 0 

h3  =  0.5  f 2. 645  =  2.145  inches. 

The  maximum  theoretical  height  of  the  diffuser  to  ob¬ 
tain  a  constant  adiabatic  Mach  number  of  2.17  could  not  be 
obtained  in  the  diffuser  section  due  to  physical  limita¬ 
tions  of  the  tunnel  itself.  As  may  be  noted  in  Figure  19b, 
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the  diffuser  contours  actually  existing  are  far  from  the 
theoretical  solution  of  the  maximum  diffuser  height*  The 
actual  contour  required  is  of  an  exponential  nature  as 
developed  previously.  Values  required  in  the  diffuser  to 
produce  a  constant  adiabatic  Mach  number  in  the  diffuser 
vs.  diffuser  length  are  presented  in  Table  XIII  in  tabular 
form  and  plotted  in  Figure  19b.  Numerical  values  of  x  were 
inserted  in  the  exponential  equation  derived  for  determin¬ 
ing  the  diffuser  height. 

i 

DETERMINATION  OF  MINIMUM  DIFFUSER  AREA  TO  ALLOW  SUPERSONIC 
FLOW. 

As  in  the  case  of  the  diverging  nozzle,  it  was  assumed 
the  diffuser  converged  linearly  from  the  maximum  nozzle  exit 
height  of  1.718  inches  to  some  minimum  height  less  than 
1.718  inches.  The  same  isoclinic  technique  was  employed  to 
determine  this  minimum  height  and  allow  supersonic  flow  to 
exist.  The  most  severe  condition  that  may  be  encountered 
is  to  have  a  normal  shock  at  the  nozzle  exit  or  at  some 
other  location  upstream  of  the  diffuser  throat*  This  con¬ 
dition  will  exist  upon  the  initial  start  of  the  tunnel  blow* 
The  minimum  diffuser  height  must  be  sufficiently  large  to 
pass  the  given  mass  rate  of  flow  in  the  subsonic  condition* 
Assuming  a  normal  shock  occurring  at  the  nozzle  exit,  with 
a  Mach  number  of  2.17  existing  in  the  nozzle  exit,  the 
initial  Mach  number  for  the  diffuser  will  be  0*551  as  taken 
from  the  Gas  Tables  (3)  for  normal  isentropic  shock.  At 
the  minimum  area,  the  Mach  number  must  increase  to  1*00* 
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Various  values  of  the  minimum  diffuser  height  h  are  assumed, 
the  isoclinic  field  drawn  for  these  conditions  and  the 
smooth  curve  faired  through  the  isoclinic  field „  Various 
values  of  h,  the  minimum  diffuser  height,  had  to  be  assumed 
to  meet  the  end  point  conditions,  namely,  the  curve  had  to 
have  an  infinite  slope  at  the  junction  of  h  and  the  Mach 
number  squared  equaling  one.  Table  XIV  presents  the  data 
for  determining  the  minimum  height  in  tabular  form  and 
Figure  2-F  is  a  plot  of  the  isdclinic  field  for  a  minimum 
diffuser  height  of  1„40  inches.,  All  isoclines  are  not 
plotted  to  reduce  the  clutter  near  the  terminal  point  of  the 
curves.  The  slope  of  dM  /dh  changes  radically  near  the  ter¬ 
minal  point  requiring  a  large  number  of  isoclines  to  be 
plotted. 

DETERMINATION  OF  MACH  NUMBER  AT  DIFFUSER  THROAT  WITH  DIF¬ 
FUSER  HEIGHT  =  1.40  INCHES  CONSIDERING  FRICTION 

Conditions  at  nozzle  throat  are  designated  by  t,  at 
the  nozzle  exit  by  2,  and  at  the  diffuser  throat  by  3. 

Utilizing  the  influence  coefficients  of  Table  8.2 
2 

multiplied  by  M  and  designating  the  product  F  as  outlined 
by  Shapiro  [2]  for  the  combined  effects  of  area  change  and 
friction 


Area  of  tunnel  =  width  x  height  =  wh 
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As  A  ^  h  for  constant  width  tunnel 


/•:  4f  dL 
f  ’  D 

hydraulic  diameter  D  = 

4  x  area 

4  wh  = 

Wetted  perimeter 

2  (w-^  h) 

5  wh 

w  f  h 

w  s  ^  inch 

g(^)  h  _  a  /? 
+h  s  +  an 


Substituting  this  value*'  of 


dM‘.FAf 


a 


Assume  that  the  diffuser  height  will  vary  linearly  in 
accordance  with  the  relationship 

h  -  my,  ■+  hQ 


where  m  is  the  slope  of  the  diffuser  and  x  is  the  distance 
down  the  tunnel. 

At 


x  =  x ^  =  4.802  inches,  h  =  h2  =  1.718  inches, 
and  at 


X  ES 

10.8  inches,  h 

=  h3  =  1 

..40  : 

inches 

m  - 

h  j  -be 

X3  " 

(0.8  - 

1/8 

4,8 

«  -  0 

0530 

dh 

-  mdx  * 

> 

d  x  ^ 

dJ/l 

m 

i 

dMz 

_  {-Em  -f - 
L  h 

4f  hk 
m 

dh 

rV  J 

dh 
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Let 


h  - 


+J]s. 


and 


p  -  UFh 


for  determination  of  F 


~h  =  'JillLi  —  /.  552  inches 

Z 


f  -  0.  0-023 


M2  -  S’.  /7  from  Figure  20 


tis-nl  =  $5  + 


4  f  FX_t±  + 


m 


(ye  +  %hi  -  hz j 


£  ^  2* 


l.55°j 


F  +  41D.00Z5),  l 
'A  -0.0530  \3.ne 


mz3  -  n2  -  -o.s2!  f  +  0.0122  rF 


Assuming  a  value  of  to  determine  M  and  the  values  of 
the  influence  coefficients  from  Table  B-6  [2] ,  a  process 
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of  iteration  may  be  made  to  find  the  correct  value  of 

__  •  «  rTt  l.(o3+  F.  H  j 

Try  A73=/.63  3  W  = — -p — —-(.TO 

Fft  -  4.TG  35 
F  -  -  /2,038 

/^7  2_  /vy  2  =  -  0t  Z2(F^  -f  <9t6>7?2/> 

3  2  » 

HZ  -  4.10  -  0.  £2!  (4.1  (=35)  -  OollZ  (/ Z.o38 ) 

3 

Mz  =  e.io 

A73  =  L  04-? 

Taj  743=  /.  0  4  )  aT  =  !•  ? 0 5 

Fa  =  4.704? 

Ff  -  -  /£./<% 

A7  2  -  4.70  -  a  ;?2>/  (4.101?)  -  O.OJ^Oe.  too) 

*5? 

M*  - 
M3  --  /.  (o4 
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In  addition  to  the  analytical  solution  for  determining 
the  Mach  number  at  the  diffuser  throat,  an  isoclinic  field 
plot  was  drawn  and  an  exit  Mach  number  obtained  by  this 
technique.  The  Mach  number  as  determined  by  the  use  of 
isoclines  was  1.71  compared  to  a  Mach  number  of  1.64  as 
determined  by  the  analytical  solution.  Table  XV  and  Figure 
3-F  present  the  information  necessary  to  plot  the  isoclinic 
field  and  a  plot  of  the  field  itself. 

DETERMINATION  OF  EXPERIMENTAL  ADIABATIC  MACH  NUMBERS  FOR 
SHOCK  FREE  CONDITIONS 

Assume  that  flow  to  the  nozzle  throat  is  isentropic  and 
that  the  flow  for  the  remainder  of  the  tunnel  system  is 
adiabatic.  From  the  Gas  Tables  for  isentropic  flow  (3)  the 
adiabatic  Mach  number,  M,  at  any  desired  point  may  be  deter¬ 
mined  for  shock  free  flow  conditions  under  the  product 


providing  A*  and  PQ  are  taken  as  the  area  and  stagnation 
pressure  when  the  Mach  number  is  unity.  P/PQ  = 

y\  p 

1  -  — —  with  pressures  PQ  and  AP  being  measured  by  the 
use  of  manometers.  Consider  pressure  tap  16  at  the  outlet 
of  the  nozzle  as  station  1  and  tap  31  in  the  diffuser  as 
station  2. 

Data  is  taken  from  Table  V  or  Figure  20  for  a  shock 
free  run  with  the  diffuser  set  at  a  constant  height  of 
1.718  inches. 
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2.  34 


(VA*)  = 

(P/PQ)1  =  0.0844 

(A/A^)  (P/Pq^  =  2.34  x  .0844  =  0.1975 

Prom  the  Gas  Tables  for  isentropic  flow  (3)  using  the  above 
value  for  entering  the  tables  M^  is  determined  as  being 
2.125.  All  adiabatic  Mach  numbers  upstream  of  a  shock  were 
determined  in  this  manner. 


DETERMINATION  OP  EXPERIMENTAL  ADIABATIC  MACH  NUMBERS  FOR 
SHOCK  IN  NOZZLE  OPERATION. 


The  adiabatic  Mach  number  for  points  downstream  of  the 

-i 

shock  in  the  nozzle  are  determined  in  the  same  manner  as 
for  shock  free  conditions.  Taking  data  from  Table  IX  for 
pressure  tap  number  20,  the  product  (A/A*)  (P/PQ)1  was  cal¬ 

culated  as  being  0.941.  Entering  Table  30  (3)  with  this 
value,  the  adiabatic  Mach  number  was  determined  as  being 
0.594.  The  ratio  of  change  in  stagnation  pressure  at  pres¬ 
sure  tap  20  or  at  any  other  location  downstream  of  the  shock 
in  the  tunnel  to  the  initial  stagnation  pressure  is  obtainec) 
by  the  relation 


fe)„ 


The  value  of  (P/P0y)n  ke;*-n9  evaluated  from  table  30  (3)  at 

the  Mach  number  Mn  under  the  column  P/PQ. 

<  , 

DETERMINATION  OF  EXPERIMENTAL  MEAN  FRICTION  COEFFICIENT 

To  eliminate  as  many  variables  as  possible,  data  with 
the  diffuser  set  at  a  constant  height  of  1.718  inches  and 
the  tunnel  running  shock  free  was  used,  i.e. ,  Fanno  Flow. 
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With  the  value  of  M  having  been  determined  from  data  in  Table 
V  as  being  2.125,  at  point  1,  the  constant  area  inlet,  this 
Mach  number  was  used  as  the  argument  to  enter  the  Gas  Tables 
(3)  for  Fanno  Flow  and  evaluate  the  following: 


(£),  -  3757 

» 

-  0.3407 

Now  since 

£  _  (P/P*L 
p  ~  (P/Pi, 


we  obtain 


( 


P 

PI 


e 


?  I0A*  (0.  3131) 
0. 844  v 


=  0.  481 


From  the  Gas  Tables  (3)  for  the  Fanno  Line  flow  using 
(p/p*) ^  -  0.481  the  Mach  number  at  the  constant  area 

pressure  tap  31,  point  2,  is  determined  as  being  1.78  and 


87 


the  corresponding  value  of 


0. 23520 


4f  a  L  \  Q'  37  37  —  0.  £3  50  -  0.1385 

O  4 

=  0.175  inches 


&  L  =  /0.0O*  -  4.75"  ^  5.25  inches 


0. 1358  »  0.775 

_5. 

£).  OJS&5 
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aoo 


K 


8  9 


h- inches  h  —  inches 

Figure  l-F.  Iaocllolc  Flald  for  Datarmloation  of  tho  Mach  Numbat  Figura  3-F.  Iaocliole  Flald  for  Datarnlnatlon  of  tho  Mach  Nt»bar 

at  Exit  of  Sagaant  Tvo  of  Mottla  Whlla  Oparatiog  Shock  Fraa.  at  tha  Hloira  Olffuaar  Araa. 
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TABLE  XI.  Simplified  Nozzle  Coordinates 
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TABLE  XIV.  Isocline  Data  for  Determination  of  Kinimum  Diffuser  Area  TABLE  XV.  Determination  of  Mach  Humber  at  Minimum  Area  Conditions  by  isoclines 
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APPENDIX  G 


AIR  RECEIVER  CAPACITANCE  DESIGN  CALCULATIONS  AND  PREDICTED 
TEMPERATURE  VARIATIONS 

To  maintain  an  isothermal  blowdown,  Muto  115]  states 

M  C7 

the  ratio  (I).  C  -  - - - ^ - -  should  have  as  a  minimum 

nb  c„ 

value  of  7.5,  and  the  non  dimensional  number  of  transfer 

h  A 


units  (2),  NTU  = 


Cv  rh 


should  have  a  minimum  value  of  7.5. 


Storage  tank  dimensions 

48  in  Outside  Diam 
11  ft.  4  in^  height 
.  75  in.  heads 
.  781  in.  shells 
Mass  in  tanks  =  Mc  Ibm. 


-  2  ?,„e*/  \jrdt/  +  ULp] 

_  r  77  -  47.Z?*  0.16/<  /it  -f 

■^WiL  me 


Mc  -  10 j  4 10  !b  mass  of  sfee! 


Surface  Area  of  tanks 


-f[-rrd.  2  R  = 

u  inside  a  J 


=?[■ 


77(4-6 .44) 
77- f  4  6.^  )/;?«&  ■+  g - 

/44 


flz 


] 


/\PEA  -  183.5  scjuar feet 


93 


Mass  of  air  initially  in  tank  as  beginning  of  blowdown  =  Mb 


Initial  tank  pressure  =  475  psia. 


£ 

Mb  =  PV  ,  iZ±  *0/3  '%*  *  g  tarts  '  m  fihanK 

RT  S3.  34  £j_jQ  x  520° R 

I  bn  -  °R 

Mb  =  577  lbm.  air  initially  in  tank. 

Mass  of  air  in  tanks  at  and  of  blowdown  at  150  psia 

Me  =  sir  =  182  lbm  @  150  P5/a 

4 15  ' 

■% 

m*  =  -  =  182/511  -  0.3 / 

^6 

Substituting  the  above  values  into  equations  1  and  2 


'  *  -  QjL  st&el  _  /Qj  4/0  X  &  U3  _  u  a 

°  ‘  Mb  C„  ~  577*  0J13  "  * 

-5  I  100  lhm/hr 

h  A  _  3  x  /<g3,.5 

Cv  0.173*1/00 

2.89 

Since  the  value  of  Cj*  7.5  and  the  value  of  NTU  ^7.5 

the  value  of  NTU  must  be  increased  by  increasing  the  area 
available  for  the  heat  transfer 


using  a  value  of  h  = 


NTU  = 


NTU  = 


NTU  = 

A  = 


7.5 


h  A 

Cv  on 


~j£_  X  O.  L23.  i  L/°°- 

3 


d7<ofis 
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additional  area  required 

476  ft*  -  /&3,5  ffZ  ■=  Z1S.5  fi* 

if  1/8  inch  aluminum  sheet  is  used  to  reduce  to  corrosion 
problem,  the  mass  of  aluminum  required  utilizing  both  sides 
of  sheet  for  heat  transfer 

Mass  A1  =  S3ZJLS£*  QJ^S  0-]  /69  [uesff3] 

Mass  A1  =  257  lbm. 

Utilizing  eqn  9A  as  developed  by  Muto  [15]  to  predict  the 
storage  tank  air  temperature  after  a  blow  down  with  the 
presently  installed  tanks  without  additional  corrective 
heating  surfaces  and  comparing  this  value  to  the  tempera¬ 
ture  actually  noted  after  a  blowdown  very  close  correla¬ 
tion  was  noted. 

Eqn  9A  by  Muto  to  predict  temp  after  blow 

*  T  I  G?)-4  -  I&3I  +  //.S) 

t0  ~  \  /+  Co  '  ~  v  /  +  t/.e  ) 

T*  =  0.9785 

with  To  =  523°R  from  Figure  1-G 
Tend  of  blow  =  0.9785  x  523°R  =  511°R 

From  Figure  1-G,  actual  temperature  of  the  storage  tank 
air  at  the  end  of  a  10  minute  run  «  510. 5°R. 
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